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Debris disks around Sun-like stars 
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ABSTRACT 



We have observed nearly 200 FGK stars at 24 and 70 microns with the Spitzer 
Space Telescope. We identify excess infrared emission, including a number of 
cases where the observed flux is more than 10 times brighter than the predicted 
photospheric flux, and interpret these signatures as evidence of debris disks in 
those systems. We combine this sample of FGK stars with similar published 
results to produce a sample of more than 350 main sequence AFGKM stars. The 
incidence of debris disks is 4.2l 2,( j ) % a t 24 microns for a sample of 213 Sun-like 
(FG) stars and 16.4+ 2 ;^% at 70 microns for 225 Sun-like (FG) stars. We find that 
the excess rates for A, F, G, and K stars are statistically indistinguishable, but 
with a suggestion of decreasing excess rate toward the later spectral types; this 
may be an age effect. The lack of strong trend among FGK stars of comparable 
ages is surprising, given the factor of 50 change in stellar luminosity across this 
spectral range. We also find that the incidence of debris disks declines very slowly 
beyond ages of 1 billion years. 

Subject headings: circumstellar matter — planetary systems: formation — in- 



frared: stars 



1. Introduction 

Planetary system formation must be studied through various indirect means due to the 
relative faintness of distant planets and the long timescales over which planetary system 
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evolution takes place. One powerful technique has advanced substantially in the era of the 
Spitzer Space Telescope: investigations of dusty debris disks around mature, main sequence 
stars. Debris disks arise from populations of planetesimals that remain from the era of planet 
formation; the analogs in our Solar System are the asteroid belt and the Kuiper Belt. Any 
system that possesses a debris disk necessarily has progressed toward forming a planetary 
system to some degree. 

The many small bodies that inhabit a debris disk can, on occasion, collide, producing 
a shower of fragments that grind each other down to dust particles. These dust grains 
can be heated by the central star to temperatures ~100 K, where they can be detected at 
wavelengths of 10-100 microns. Dat a from the IRAS and I SO satellites were used to identify 
and characterize debris disks 



200ll ; lHabing et al. 



2001 



(e.g., 



Aumann et al.l Il984j ; iDecin et al.l 120001 ; iSpangler et al. 



Decin et al.ll2003l ). The Multiband Imaging Photometer for Spitzer 



(MIPS; iRieke et al.l (l2004j )) offers substantially improved sensitivities at 24 and 70 microns 
and therefore can be used to advance the study of debris disks and measure the fraction of 
stars that possess colliding swarms of remnant planetesimals. 



A number of these surveys have been carried out using MIPS. IRieke et al.l (120051 ) and 



Su et al.l (120061 ) observed hundreds of A stars and found that the number of A stars showing 
thermal infrared excess suggestive of collisionally-produced dust decreases as a function of 
stellar age, from 50% or more at ages just past the gas dissipation age of 10 Myr to 3 0% at 
500 Myr. The excess rates are lower for older and lower mass stars. iBryden et al.l (120061 ) 
found that the excess rat e is around 1 5 % for stars of mass and age similar to our Sun 
(for a sample of 69 stars). iGautier et al.l (120071 ) did not detect any 24 or 70 micron excesses 
suggestive of debris disks in a sample of 60 field (old) M stars (with only 13 strongly detected 
at 70 mi crons). Finally, mul tiplicity appears to play a significant role in modulating debris 
disks, as iTrilling et al.l (120071 ) found that the debris disk rate for A and F binaries was higher 
than that for single stars of similar spectral types. 

We present here a survey for excesses around almost 200 F, G, and K st ars, with ages 



and m asses similar to that of our Sun. (Some of these data were published in IBryden et al. 



( 120061 ).) We characterize the excesses that we find and present a few systems of particular 
interest. We create a larger sample of FGK stars by adding 75 stars from a similar survey, 
and derive excess rates as a function of spectral type and as a function of age. Finally, we 
discuss the implications of our results for planetary system formation. 
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2. Sample selection 

The data presented here originate in two separate observing programs. The majority 
are from a survey targeting nearby late F, G, and early K "Solar type" stars (Spitze r GTO 



program #41). Results for 69 stars from this program were published in lBryden et al.l (120061 ) 
(hereafter B06), and the present paper completes the analysis of data from that program. 
For completeness, all targets that appeared in B06 are also included in this paper. The 
remainder of the data belong to a survey for debris disks among F stars (Spitzer GTO pro- 
gram #30211). Stars in this spectral type range were omitted from previous Spitzer disk 
surveys, and have particular relevance in determining the nature of the transition from the 
high A star debris disk incidence to the more modest incidence for Sun-like stars. F urther- 



more, these 49 F stars will serve as a useful control sample to the iTrilling et al.l (120071 ) result 
that binary A and F stars have a relatively high excess rate. We combine the data from 
these two samples here to give continuous coverage, a large sample of stars, and overall view 
of debris disks from F0 to K5 stars. We refer to the combined sample of 193 stars (PID 41 
plus PID 30211) as our FGK sample. Relevant information for the combined target list is 
given in Table [TJ Figures dHU show the distributions of targets as functions of spectral type, 
distance, age, and metallicity. 



2.1. F stars sample 

Our target list was assembled in the following way. A list of main sequence F0-F5 stars 
was created from the Hipparcos database, in order of increasing distance from the Earth. 
All stars that were known to be members of multiple systems were removed (by vetting 
against the CCDM catalog). Next, system s that were obvio u sly yo ung (less than 1 Gyr) 



were removed, using the measurements of iNordstrom et al.l (120041 ) . The average age of 



our sample is around 3 Gyr (compare to the average age of the "Solar type" FGK sample 
presented in Section [2T2l of around 5 Gyr). We used the IPAC IBIS tool to extract the infrared 
background and confusion levels from the IRAS data for each remaining target, compared the 
predicted target photospheric flux to the confusion per MIPS beam, and discarded targets 
with predicted S/N<3. Finally, we discarded targets requiring 70 micron integration times 
longer than 1300 second (12 cycles of 10 second integrations); targets fainter than this rapidly 
become very expensive with questionable return. These selection criteria left 49 targets as 
our final target list. Here we present the 37 of those targets observed to date; the remaining 
stragglers will be published in a future paper. 
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2.2. "Solar type" sample 

The "Solar type" sample selection is described in detail in B06. We briefly summarize 
that description here. This program consists of two overlapping sets of stars: those that 
meet strict selection criteria for an unbiased sample, and those that are known to harbor 
planets. In both cases, only stars with spectral type similar to the Sun are considered. Stars 
with spectral type F5 to K5 and luminosity class IV or V were included, with the final list 
chosen by expected signal to noise ratio for photospheric flux. Additionally, a minimum 
photospheric 70 micron flux was set for each spectral type bin: 20 mJy for F5-F9 stars, 10 
mJy for G0-G4, and 5 mJy for G5-K5. There is no explicit selection based on stellar age or 
metallicity; however, the 70 micron brightness and S/N thresholds are relaxed in some cases 
to allow stars with well determined ages into the sample. Close, resolved binary systems are 
excluded. The total sample pres ented here is 156 s tars (non-planet stars and planet stars, 



combined). A companion paper (IBryden et al.l 120071 ) describes the planet-bearing sample in 
more detail. 



Observations and data reduction 



A listing of the observations for these programs is given in Table [2j All Spitzer obser- 
vations reported here were made between January, 2004, and March, 2007. We used MIPS 
to observe each system at 24 fim and 70 /j,m (effective wavelengths 23.68 and 71.42 /im, 
respectively). All stars were observed in MIPS Photometry small-field mode. The 24 /im 
observations all used 3 sec DCEs (data collection events) and a single template cycle (with 
a few exceptions; see Table |2]). The 70 /im observations typically used 10 sec DCEs and 5 to 
10 template cycles. 



Ove rall, our data analy sis is similar to that previous ly described in iBeichman et al. 



( ]2005al ) , IBryden et al.l (120061 ) , and IBeichman et al.l (l2006bl ) . The data reduction is based 
on the DAT software developed by the MIPS instrument team ( Gordon et ali 2005 ). At 
24 mic rons, image flats are chos en as a function of scan mirror position to correct for dust 
spots (jEngelbracht et al.ll2007al ). Images were mosaicked from individual frames with half- 
pixel subsampling. We carried out aperture photometry on reduced images as described in 
Beichman et al.l (J2005aJ), using an aperture of six 2/55 pixels, a sky annulus of 12-17 pixels, 
a calibration factor of 1.047 /iJy/arcsec 2 /(DN/s), and an aperture correction of 1.15. At 
70 microns we used images processed beyond the standard DAT software to correct for 
time-dependent transients, correcti ons that can significantly improve the sensitivity of the 
measurements (IGordon et al.ll2007l ). Because the accuracy of the 70 micron data is limited 
by background noise, rather than instrumental effects, a very small photometric aperture was 
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used to maximize signal-to-noise - just 1.5 pixels in radius, requiring an aperture correction 
of 1.79. An annulus with an inner radius of 4 pixels and an outer radius of 8 pixels was 
used to measure the sky background around each target. The flux level is calibrated at 
16.5 mJy/arcsec 2 /MIPS_70_unit, with a default color correction of 1.00 (MIPS_70_units are 
based on the ratio of the measured signal to the stimulator flash signal). For both the 24 and 
70 micron data, neighboring point sources were subtracted from the images before measuring 
the sky brightness. 

The 24 mic ron centroid posit ions, which are consistent with the telescope pointing 



accuracy of <1" ([Werner et al.ll2004j ). are used as the target coordinates for both wavelengths. 
The only significant update from previous papers is a 4% increase in the overall calibration 
at 70 microns from 15.8 to 16.5 mJy/arcsec 2 /MIPS_70_unit. 



4. Photospheric predictions and detections of excess 

4.1. 24 /xm 

To determine whether any of our target stars have an IR excess, we compare the mea- 
sured photometry (F) to predicted photospheric levels (P). High quality K magnitudes are 
required to extrapola te to 24 microns. Most of our stars have good 2MASS photometry 
( iSkrutskie et al.ll2006l ). but stars brighter than mx ~ 4.5 saturate in 2MASS and therefore 



this survey cannot provide accurate photometry. We obtained photometry for as many of 



Bessell et al. 


(1998 


) and 


Carpenter 


(2001) 



some of this "heritage photometry" are difficult to estimate, since infrared photometric sys- 
tems were still under evolution when some of it was obtained. We show below that the net 
errors after transformation are surprisingly small. For 16 additional stars where we could 
find no K-band photometry, we estimated values from B-I or B-V colors. We obtained the 
optical colors from SIMBAD for the members of our sample and carried out linear fits against 
I-K or V-K. The rms scatter around the fit was 0.07 for B-I vs. I-K and 0.09 for B-V vs. 
V-K, and we take this scatter to be the error in the estimated K and K5 magnitudes. There 
are 12 stars for which we were unable to make even estimates of K and Kg magnitudes. 

For the 181 sources where K5 data is available, we use the following approach. We 
computed ratios of the measured flux densities at 24 microns and K5 to determine the mean 
photospheric color. There is a weak dependence of this photospheric color on spectral type. 
We calculated this correction by assigning a numerical value proportional to spectral type (0 
for AO, 10 for F0, 15 for F5, etc.) and multiplying the ratios by 1.048/(l+0.048*(type/22)), 
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determined from a linear fit to the data. This procedure results in a predicted 24 micron 
flux (P24) for each of these 181 stars. 

We then calculate the excess ratios (R24=F24/P24) for each system (Table [3D- The 
distribution of derived excess ratios is shown in Figure [51 The R24 distribution for all stars 
with Ks magnitudes (2MASS or heritage) is quite regular, with clear peak near unity. For 
0.92<R24<1.08, the mean and median both have a value of 0.995, with an rms scatter 
of 0.032. For only the stars with heritage photometry, the analogous average is 0.997, the 
rms scatter is 0.031, and the median is 0.999. That is, there is no evidence for any systematic 
difference in the sources of K photometry and we can treat all stars measured in this way 
identically. We can therefo re set a 3a thres hold for identification of an excess. Our approach 



is the same that we took in lSu et al.l (120061 ). We take la to be the rms scatter value of 0.032. 
Our result of zero sources with R24<0.92 is consistent with this value of a. We therefore set 
the 3a excess threshold to be 1.10. One star (HD 222404) has R24=1.10, and we consider 
this marginal case explicitly in the next paragraph. The next smallest R24 value above our 
3a threshold is 1.15 (close to 5a), and six stars have even larger R24 values. The stars with 
2MASS Ks magnitudes or heritage K magnitudes and excesses at 24/zm are HD 166, 3126, 
69830, and 105912. 



The stars with possible excesses based on estimated K magnitudes need to be discussed 
individually because of their larger K magnitude uncertainties. HD 10647 has a ratio of 
1.21, based on extrapolation from B-V and hence at a level of ~2.5o\ It has a very large 



excess at 70 /im, by more than a factor of ten over the photosphere, and I Chen et al.l (J2006) 
find an excess at 30-34 microns, so we consider the excess at 24 /zm to be confirmed at the 
longer wavelengths. HD 101259 has an indicated excess ratio of 1.25, based on extrapolation 
from B-V, and hence at nearly 3a. It is not saturated in the 2MASS data at J, and if we 
apply the standard J-K color for its spectral type, we estimate an excess ratio of 1.15. It 
has a modest but significant excess at 70 /im, so we accept the reality of the indicated one 
at 24 (Ltm. HD 191408 has a ratio of 1.23, based on an estimate from B-I and hence at a 
level of ~3cr. It has no significant excess at 70 /im, with a 2a upper limit to the excess 
ratio of 1.7. We used Kurucz model fitting to make an independent prediction of the 24 /im 
flux density, and found an excess ratio of only 1.08, suggesting that our extrapolation is in 
error, probably due to a bad photo metric point. Furtherm ore, IRS spectra throughout this 



region are photospheric (no excess) (IBeichman et al.ll2006al ). It is therefore unlikely that the 



possible 24 /im excess is real. HD 222404 has an excess ratio of 1.10 based on a K magnitude 
estimated from B-I, and hence at a level of 1.5er. However, the 2a upper limit to the excess 
ratio at 70 /im is only 1.04, so we reject the possibility of an excess at 24 /im. The remaining 
stars with estimated K magnitudes all have excess ratios <1.10. In summary, we therefore 
find six stars in the sample of stars with K magnitudes (2MASS, heritage, or estimated) 
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with significant 24 /im excesses (Table EJ: HD 166, 3126, 10647, 69830, 101259, and 105912. 

This "self calibration" technique (using the data to define the outliers) is generally re- 
quired in determining intrinsic stellar colors. For extreme outliers (Figure ED, there is no 
question of the veracity of this technique. Furthermore, because the distribution of R24 is 
well-behaved, with very small scatter, we are confident that it robustly reflects the true dis- 
tribution of excesses. We further note that there is no procedure other than self calibration 
that works. There is no sample of stars measured at 24 /jm that includes the relevant spectral 
types and that is certain to have no sm all excesses. Ku rucz model fitting leaves a small resid- 
ual that must be removed empirically (ISu et al.l 120061 ) . We checked our results using model 
fitting for a subset of targets and found agreement with the photometry-based extrapolations 
performed for the entire sample. Finally, we note that the rms scatter in the distribution 
of R24 is 0.02-0023 from the K-band photometry, 0.01 from the 24 jum photometry, and 
about 0.01 from the color-color relationships. Assuming these errors add quadratically (they 
are independent errors), the error left for variations in stellar color is 0.019, a very small 
number. This leaves little room for significant excesses for the bulk of the stars. We conclude 
that our self calibration technique is robust. 



4.2. 70 /im 

Since the majority of the stars in our sample do not have 24 micron excesses, in general 
we make our 70 micron photospheric predictions based on the 24 micron flux measurements. 
For stars with no excess emission, the expected F24/F70 ratio is 9.10 (scaling by A 2 and 
using the effective wavelengths specified above), so we derive predicted 70 micron fluxes 
by dividing the observed 24 micron fluxes by this numerical factor. The predicted fluxes 
(P70) are listed in Table [3j Calculating P70 this way (and not from extrapolations from 
K magnitudes) allows us to capitalize on the nearly uniform 24/70 flux ratio rather than 
extrapolating from K and potentially magnifying any slight misestimates. It also allows us 
to determine P70 for the 12 stars that have no K band measurements and therefore no P24, 
so that we can search all 193 stars in our sample for 70 micron excesses. Naturally, for the six 
stars with 24 micron excess this extrapolation from observed 24 micron fluxes to 70 micron 
predictions does not work; for these systems, we extrapolate the predicted 24 micron fluxes 
to derive predicted 70 micron fluxes, using this same numerical factor of 9.10. For these six 
cases, our 70 micron predictions are confirmed with Kurucz model fitting. 

The large scatter and large number of (expected) excesses at 70 microns makes deter- 
mining the excess threshold in the circular manner above problematic. We instead quantify 
the significance of the (possible) 70 micron excess X70> which is defined as follows: 
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F70 - P70 

X70 = 

O"70 

where F70 and P70 are as defined above and aro is the error in the measurement. 



We use the same criterion as ISu et al.l (120061 ) to identify significant excesses: we require 



X70 > 3. This is identical to saying that the significance of the excess must be at least 
3a. The distribution of X70 is shown in Figure [6j The core of this distribution is clearly 
centered near zero, the expected value for a sample where non-excess is the majority outcome. 
Nevertheless, the distribution of X70 shown in the top panel of Figure E] is not symmetric in 
the range —3 < X70 < 3: the broad shape of the right "shoulder" in the top panel of Figure E] 
(1.5 < X70 < 3) suggests the presence of a number of low-level excesses (below our detection 
threshold) . 

In our large sample, we can test the reliability of these determinations by counting the 
number of stars that fall below our -3a excess criteria (that is, X70 < —3.0). Only one star 
(HD 100067) has X70 < —3.0 (Figure [6]). This system has a large negative xm value because 
the measured aperture photometry is formally negative (meaning the sky value subtracts 
more than the enclosed flux in the photometric aperture); in other words, this star was 
simply not detected at 70 microns (Table ED- There are no cases of X70 < —3.0 where the 
star is detected. Assuming symmetric noise properties, we therefore expect no spurious 
positive outliers. 

In summary, we identify any system with X70 — 3.0 as an excess system. There are 
27 such systems out of the 181 stars with K magnitudes. We define HD 69830 to have an 
excess at 70 microns because it has X70 = 2.81 and a significant (corroborating) excess at 
24 microns (R24=1.47), bringing the number of systems with 70 micron excesses to 28 out 
of 181. This 28 includes 6 systems that have excesses at both 24 and 70 microns. Finally, 2 of 
the 12 stars with no K magnitudes also have 70 micron excesses (see Section |5\3"I) . based on 
F24/F70<9.10. This brings the total number of systems with 70 micron excesses to 30 out 
of 193 (Table [3J . 

In Figure [7] we show the distribution of our 181 sources with K magnitudes in R24-X70 
space. We confirm that there are likely no other spurious 24 micron excesses. Furthermore, 
the random scatter distribution (fully uncorrelated) of the non-excess sources confirms that 
there are no systematic errors introduced by our technique of predicting 24 micron fluxes 
and extrapolating to predict 70 micron fluxes. Systems with weak excesses at both 24 and 
70 microns would not necessarily be identified by the techniques we have described here. 
However, we see in Figure [7J that there are no systems in the upper right of the enclosed 
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area, where such a system with two weak excesses would reside. 



Results 



5.1. Excess rates 



There are 181 stars in our sample with well-constrained K band magnitudes and there- 
fore reliable 24 micron predicted fluxes. We can make excess determinations for all 193 stars 
at 70 microns (in most cases based on our 24 micron measurements). We add to this catalog 
three stars that meet the selection criteria for PID 41 but had already been included in 
other observing programs: e Eri (K2V, 3 x 10 8 yr ); r Ceti (G8V, 7.2 x 10 9 yr); and r 1 E ri 
(F6V, 3 x 10 s yr) (Table 3 from B06 as well as becin et al.l pS); Ichen et al.l (bood )). 
The first of these has excesses at both 24 and 70 microns (Backman et al., in prep.); the 



second is photospheric at 24 microns and has an excess at 70 microns (IChen et al. 



2006 



Stapelfe ldt et al.ll2006l); and the third has no ex cess at either 24 or 70 microns (IChen et al. 



2006 



Stapelfeldt et al.ll2006uTrilling et aU2007l ). These additions bring the total excess rates 



to the following values: at 24 microns, 7 out of 184 stars have excesses, giving 3.8_{^%; and 
at 70 microns, 32 out of 196 stars have excesses, giving 16.3l|g%. In all cases in this paper, 
we cite binomial errors that in clude 68% of the proba bility (equivalent to the la range for 
gaussian errors) , as defined in iBurgasser et al.l (120031 ) . All of the systems with 24 micron 
excesses also have 70 micron excesses. These excess rates are also presented in Table HI 

Of the 196 stars presented here in our FGK sample, 48 are known to harbor extraso- 
lar planets. The excess rates for the planet sample are higher than the no-planet sample, 
although formally t he two samples' excess rates are consistent at the la level (Table H|). 
Bryden et al.l (120071 ) discuss further the differences for excess rates and debris disks for stars 
with and without known planets. Planets have be en found around ~10% of the stars sur- 
veyed in radial velocity programs (jUdry et al.ll2007l ). but here planet-bearing stars represent 
25% of the total FGK sample, so there is the potential for contamination of our derived 
excess rate for FGK stars. However, if we weight the planets and no-planets excess rates 
(Table UJ) by 10% and 90%, respectively, we find that the resulting total excess rates are 
3.9% (24 microns) and 15.6% (70 microns), hardly different from the overall excess rates 
derived above (3. 8^2% and 16.3+2!%, respectively). We conclude that any contamination 
must therefore be insignificant. 
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5.2. Dust properties 

We interpret the presence of excess thermal infrared flux as a signature of emission from 
dust grains. These dust grains are assumed to be produced relatively recently from collisions 
of asteroid-like bodies. (Dust grains in these systems have short lifetimes against radiation 
forces.) By constraining the properties of this dust, we can learn about the processes of 
planetary system formation in a large sample of FGK stars. 

In general, we would like to fit the observed excesses to a black body to determine the 
temperature of the excess. For systems with excesses at both 24 and 70 microns this is rel- 
atively straightforward: we find the blackbody temperature that best fits the two observed 
excesses. These temperatures are given in Table [5j We scale this blackbody by the appropri- 
ate amount to fit the excess measurements. The ratio of this scaled blackbody's integrated 
flux to the star's integrated flux is the fractional luminosity (Table Figure [8]). 

Most stars in our sample that show excesses, however, have excesses only at 70 microns, 
which leaves dust temperatures relatively unconstrained. In these cases — excess flux at 
70 microns, photospheric flux at 24 microns — we set the 24 micron "excess" to be equal to 
three times the error at 24 microns (which is dominated by the calibration uncertainty), and 
find the blackbody that best fits this "excess" and the measured 70 microns excess. This 
resulting temperature is the maximum temperature for the excess (Tabled]). 

After solving for this dust temperatu re (2V), we can calcu l ate th e orbital distance r (in 



AU) of the dust through equation 3 from lBackman &: Parescd (Il993f ) : 



278 \ 2 f U x 



Td J \L e 



(1) 



where is the stellar luminosity. (Large grains with unit emissivities that radiate as black 
bodies are assumed.) We calculate the stellar luminosity of the host star simply through 

U = 4tt2^ / Kurucz model = AnRloT^ ff (2) 



where R* is the stellar radius (from iDrilling fc Landoltll2000l ) and the stellar effective tern 



perature T e ff is given in Table [T] and again in Table El Our calculated dust distances are 
given in Tables [5] and [6] and Figure El For systems with excesses only at 70 microns, the dis- 
tances we derive are minimum distances, since the dust temperatures we derive are maximum 
temperatures. In all cases, we assume black body grains. 

Finally, from our blackbody fits to the excesses we also derive the fractional luminosities: 
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the ratios of dust luminosity to stellar luminosity. This quantity is also referred to as fd and 
Ld/L*. The fractional luminosities we derive are given in Tables [5] and [6] and Figure [8] As 
above, for systems with excesses only at 70 microns, these derived fractional luminosities are 
the maximum values that are consistent with the observed data. However, extremely massive 
cold outer disks, which could imply a larger fractional luminosity than the "maximum" values 
we present here, cannot be ruled out. 

To place a reasonable lower limit on the fractional luminosity e stimates for systems 
with excesses only at 70 microns, we follow the arguments in B06 and iTrilling et al.l (120071 ) 
and assume a "minimum" possible disk temperature of 50 K. This assumption excludes 
the possibility that a very cold (and potentially massive) disk could exist, with emission at 
wavelengths longer than 70 microns. This 50 K blackbody is scaled to match the measured 
70 micron excesses (this results in no detectable excess emission at 24 microns, consistent with 
the observations). We calculate the ratio of this disk's luminosity to the stellar luminosity, 
to derive the minimum fractional luminosities for systems with excesses at 70 microns only. 
We present these minimum disk luminosities in Table [6] and Figure [8j These two estimates 
of fractional luminosity give the range of values that fit the data. In many cases, the range 
of acceptable values is quite small. We see that th e minimum fractional luminosity is arou nd 
10~ 5 , as has been found in other surveys as well (IBryden et al.ll2006t ITrilling et al.l 120071 ) . 



5.3. Individual systems of interest 

We describe here a few interesting individual cases, some of which are also shown in 
Figure [HI to demonstrate the range of disk properties evident in our sample. 

HD 693. This system has Xm — 2.95, just below our excess threshold. It is likely that 
there is a real excess in this system that is simply slightly too faint for us to detect formally. 
Additionally, a massive and/or cold(er) disk cannot be ruled out. We identify this system 
as one of particular interest because of its low fractional luminosity: 3 x 10~ 6 (assuming 
that there is a real excess), which would make it comparable to some estimates of our Solar 
System's Kuiper Belt (Figure E]). Also of note is its dust distance of 3.2 AU, clearly well 
within the realm of possible planetary systems. 

HD 3126. This star has a very large 70 micron excess (X70 = 35.8) and a modest 24 micron 
excess (R24=1.16). With an age of 3.5 Gyr and a SIMBAD spectral type of F2, there is 
nothing apparently anomalous about this system that co uld explain its large excess, i.e., the 



star is not young, nor is it evolved. (IFavata et al.l (119931 ) assign a spectral type of F4V.) In 



Figure [9], we show the narrow range of possible fractional luminosities (~10 4 ) that are in 
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accord with the observed data. The fractional luminosity is relatively high. This system is 
clearly promising for additional observations to better characterize the dust population. 

HD 10647. This star has a well-known infrared excess that was first detected with IRAS 



(jStencel &: Backmanlll99ll ). I Chen et al.l (120061 ) determined a fractional luminosity of 15 x 
10~ 5 for this system with Spitzer/IRS measurements, using a temperature of 70 K. This 
is fairly consistent with the values we derive (Table ED, and is a relatively large fractional 
luminosity for the FGK sample we present here. HD 10647 is also known to harbor (at least) 



one extrasolar planet, orbiting at 2 AU (IBut 



er et al.l 



Jura et al. 



2006). We find the dust to be at a 



(120041 ) find an inner disk radius of 



distance of at least 21 AU (Table ED, although 
11 AU using IRS spectra and a warmer excess temperature than we find. This disk has also 
been detected in scattered light (Stapelfeldt et al., in prep.). This system is clearly fertile 
ground for further studies of planetary system formation, since both a giant planet and a 
large dust disk (and implied collisions among asteroidal bodies) are known to exist. 

HD 19994. There is no good K magnitude for this system, so we cannot determine whether 
it has a 24 micron excess. However, F24/F70 is 5.6, making this system far too bright at 
70 microns for the [24:70] color to be photospheric; we determine X70 — 4.25. This assumes 
that F24 is photospheric, which is the conservative assumption that makes our calculated X70 
a lower limit. Under this assumption, the maximum dust temperature is 170 K, the minimum 
dust distance is 4.1 AU, and the maximum fractional luminosity is 4 x 10~ 5 . These properties 
are similar to the other systems with excesses at 70 microns only. 

HD 30495. Among systems with no 24 micron excesses, HD 30495 has the most significant 
70 micron excess (%7o = 26.9). These parameters suggest a relatively cool dust disk, as seen 
in Figure O As for HD 3126, a massive, cold disk cannot be ruled out. 



HD 69830. This system was discov ered by iBeichman et al.l (l2005bl ) to have substantial 
hot emission, which iLisse et al.l (120071 ) attributed to the recent breakup of a large asteroid 
in that system. The excess infrared emission, evident from 8-35 microns (and, as we show 
here, marginally evident in photometrjQ at 70 microns), suggests quite high temperatures 
(~400 K) and a dust distance around 1 AU. This system further became interesting with 
the subsequ ent discovery of t hree Neptune mass planets in that system, all orbiting at less 
than 1 AU Lovis et al.ll2006h . 



The 70 micron flux reported here [26 mJy] is higher than that in IBeichman et "all d20Q5bh [19 mJy] 
primarily due to better centering of the photometric aperture on the source. There is an additional small 
[5%] increase in the calibration factor that has been applied since IBeichman et al] (|2005bl ). 



13 



HD 82943. Among systems with excesses at 70 microns only , HP 82943 has the largest 
fractional luminosity, at more than 10" 4 (Table E]) . Moor et al.l (120061 ) find that all systems 



with debris disk fractional luminosities greater than 5 x 10 -4 are young. HD 82943 has an 
age of more than 4 Gyr (Table [I]). Our fractional luminosity estimates for the disk around 
HD 82943 are below this apparent limit, so the age of more than 4 Gyr (Tabled]) offers no 
direct contradiction. Nevertheless, this system may be u seful in helping to define the exact 



border of the distribution discussed by lMoor et al.l (120061 ) . (Other high fractional luminosity 



systems given in Tables E] and E] may also provide useful constraints.) We note that HD 82943 



Mayor et al. 


2004; 


Lee et al. 


2006) 



HD 101259. We find an excess at both 24 and 70 microns for this star. Unusually, the 
excess in this system is not much brighter, relatively, at 70 microns than at 24 microns 
(R24=1.25 and R70=1.61). However, we note that both P24 and P70 for this system are 
based on the B-V color, and therefore have a large uncertainty. Taken together, a small and 
warm excess is likely, but not certain. Follow-up observations will be critical in characterizing 
this system, and will be of particular interest since systems with warm dust are rare. 

The measured excesses of R24=1.25 and R70=1.61 require a fairly warm disk temper- 
ature of 271 K, giving a dust distance of around 1 AU (Table E]). (The IRAS 12 micron flux 
also appears to be ~25% higher than the expected photospheric flux.) The properties of this 
disk appear quite similar to HD 69830 (see Figure [7]), and further observations (including IRS 
spectroscopy) should be made to characterize the dust population. In temperature and dis- 
tance, the dust in this system closely resembles the zodiacal dust in our Solar System, which 
is produced largely by collisions in the asteroid belt. The architecture of that planetesimal 
(and potential planet) system may not be that different from our own Solar System's. 



HD 207129. There is no good K magnitude measurement for this star, so we cannot 
predict its 24 or 70 micron fluxes. N evertheless, th is star has been known to have an in- 
frared excess for more than 20 y ears (lAumannl 119851 ) . and has been studied with both ISO 
(jJourdain de Muizon et al.lll999l ) and in a number of ongoing Spitzer studies. The F24/F70 
ratio of 0.57 that we observe here is far from the photospheric value of 9.10 and clearly indi- 
cates a strong 70 micron excess and presence of a debris disk. We fit published visible and 
near-infrared data with a Kurucz model to predict the 24 and 70 micron fluxes (Figure [TOl) 
and constrain the dust properties. 



We find P24=139 mJy and P70=15 mJy This gives excesses at 24 and 70 microns of 
25 and 274 mJy, respectively. The fractional luminosity is around 10~ 4 , which is very high 
for a star like our Sun. The color temperature is 72 K, which, for this GO star, gives Rdust of 
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15.3 AU. The temperature we derive is hotter than t he ISO-derived range of 15-45 K , but 
our fractional luminosity agrees well with their values (jJourdain de Muizon et al.lll999l ). The 
properties we derive are similar to the other systems with excesses at both bands (Tabled]). 
This debris disk ring has also been detected directly in scattered light (Krist et a!., in prep.). 



6. Discussion 

6.1. Metallicity and excesses 

Four stars in our FGK sample do not have published metallicities. For the 189 stars 
with known metallicities in Tabled], the mean metallicity is -0.08±0.22, with a median of - 
0.06. The mean metallicity of stars with excesses is -0.11±0.19 (median is -0.09). The 
mean metallicity of stars with no excesses is -0.08±0.22 (median is -0.05). There is no 
difference between the metallicity of the population of stars with excesses and the population 
of stars with no excesse s. This lack of correlation has been discussed previously (e.g., BOG, 



Beichman et all (I2006bl )). 



6.2. Excess rates across spectral type 



Many parameters affecting infrared excesses change with spectral type, including the 
importance of grain loss mechanisms (winds, Poynting- Robertson drag, photon pressure); 
stellar luminosity; and the locations of key temperatures (e.g., the ice line) in the systems. 
Significant surveys for debris disks across spectral types A-M have now been published, and 
we use those data to look for systematic trends across spectral types. 



I t is well known that excess rates decrease with stellar age (IHabing et all200ll ; ISpangler et al. 



2001uRieke et al.ll2005l ; ISu et al.ll2006l ; ISiegler et al.ll2007l ). This dependence must be avoided 
in testing for c hanges with spectral type. We take the oldest (>600 Myr) A stars from the 



Su et al 



(120061 ) sample as our representative sample from that spectral typ e. For our F, G, 



and K samples we take the union of the data presented h ere and the data i nlBeichman et al. 
( 2006b ). The targets and data reduction presented in Beichman et al. ( 2006bl ) are quite 



similar to the selections and techniques we have employed here, which allows us to merge 
the two samples relatively seamlessly. We calcul ate the excess r atios for these F, G, and K 
samples and take the M stars excess rates from iGautier et al.l (120071 ). This compilation is 
presented in Table H] and Figure [TT1 



Within the error bars, the excess rates for the A, F, G, and K subsamples are essentially 
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indistinguishable. However, there is a suggestion of a trend of decreasing 70 fxm excess 
rates with later spectral types (Figure [TTj) . We note, however, that the mean age for the 
populations increases with later spectral types. It is possible th at we are instead d e tecting a 



time-r elated effect, althou g h the decay timescales identified by iRieke et al.l (120051 ) , ISu et al. 



(120061 ) . and ISiegler et al.l (120071 ) of hundreds of millions of years should long since have 
diminished all disks at ages of billions of years. We discuss this possibility in Section | 



Beichman et al.l (l2006bl ) remarked that the excess rates for K stars appeared to be 
lower than that for F and G stars, finding zero excesses among 23 stars later than K2 
(and zero excesses amo ng a larger combined sample of 61 K1-M6 stars). We include the 
Beichman et al.l (l2006bl ) data in our analysis here, and find that, formally, the excess rate 
for K stars is no t signifi cantly different than the excess rates for earlier (F and G) stars. As 
Beichman et al.l (l2006bl ) note, and we confirm, none of the 6 K stars with excesses in our 
larger sample are later than K2. 

Part of the motivation for assembling the F stars program (PIP 30211) was as a control 
sample for the binary star program presented in lTrilling et al.l (120071 ). In that study, 69 A3- 
F8 binary star systems were found, overall, to have relatively high excess rates: 9% at 

24 microns and 40% at 70 microns. It is clear from our results here (see Figure [TTj) that 
the (single) F stars excess rate is equal to or lower than the (single) A stars excess rate. 
The binaries excess rates remain significantly high compared to the control sample of A and 
F stars. 

Figure [8] shows that, in general, there is no trend of dust distance as a function of 
stellar effective temperature (spectral type). This may suggest that the processes that drive 
planetesimal formation do not depend strongly on a single critical temperature, as would be 
the case in the "ice line" model, where protoplanetary disk surface densities increase across 
certain temperature boundaries. However, our method for calculating dust distances may 
be too crude to see this effect. We note that all of the (minimum) dust distances given in 
Table [6] would fall within the planetary realm of our Solar System (<30 AU). We are not 
observing disks that are far outside of the potential planetary realm of these systems. 

We stated above that there is no particular trend for either fractional luminosity or 
dust distance with spectral type, but there is an important caveat: no disks with large dust 
distances or relatively small fractional luminosities were identified among the latest stars in 
our sample. This dearth may simply allude to the fact that later stars are cooler. Dust at 

25 AU around a Kl (5000 K) star would have a temperature around 50 K; this dust would 
have its peak emission near 70 microns, and cooler (more distant) dust would have its peak 
emission longward. MIPS 70 micron observations of such a dust population, or a cooler one, 
would not readily show the presence of this excess. The lack of distant disks around K stars 
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may therefore be an observational bias. Similarly, low fractional luminosity disks would be 
more difficult to detect around K stars than around earlier stars, so the lack of low fractional 
luminosity disks for later stars may also be due to observational bias. These observational 
biases may be corrected with sufficiently sensitive measurements at ~100 microns (e.g., with 
Herschel) . 

To further address the question of whether there is any detectable trend of excess rate 
as a function of spectral type using existing published data, we compared the incidence of 
70 /jm excesses across these 5 spectral types (A, F, G, K, M), a total sample size of more 
than 350 stars. To avoid uncontrolled selection effects, we confined the comparison to stars 
that would have been detected at 70 microns at a level of at least 2:1 on the photosphere. We 
then applied a number of tests. First, we computed weighted average values of R70 (observed 
flux over predicted flux) as in iGautier et al.l (120071 ). We find a value of ~5 for the A stars, 
but we discount this large average because of the small size of the sample (27 stars). The 
values for the F, G, K, and M stars are 1.16, 1.23, 1.06, and 1.025, respectively. Errors are 
difficult to estimate because the excesses are not normally distributed. We also computed 
straight averages of R70 for these same samples. Here, we calculate 4, 2.6, 1.8, 1.4, and 1.1 
for the A, F, G, K, and M stars, respectively. Again, the values need to be interpreted with 
caution because error estimation is difficult. To circumvent the difficulties in determining 
errors, we binned the excesses into intervals of 0.5 in excess ratio and used the K-S test to 
determine if the resulting distributions were likely to have been drawn from the identical 
parent distribution. For each stellar type, we tested the relevant distribution against the 
distribution for all the stars, excluding the contribution of the spectral type in question. 
The result was a set of probabilities of 0.03, 0.8, 0.3, 0.3, and 0.01 that the types A, F, G, 
K, and M, respectively, were drawn from the same parent distribution as the other types. 
(For this test, a claim of a significant difference requires a probability of 0.05 or less that 
the samples are from the same distribution. Values of 0.3 imply that the the G and K stars 
are la different from their respective control samples.) From this suite of tests, we conclude 
that the incidence of excesses is different for old A stars and for M stars from that of the rest 
of our sample. We further deduce that there is a possibility of a difference appearing in the 
K stars from their lower average excess ratio. The distributions for F and G stars appear 
to be indistinguishable with our data. We employ this conclusion in the creation of an FG 
"supersample" (Section 16.41) . 

The higher excesses indicated for the old A stars could be an age effect, since the sample 
is by necessity significantly younger than t h e late r typ es (which we select ed in general to 
be > 1 Gyr in age). In fact, iGorlova et al.l (120061 ) and ISiegler et al.l (120071 ) compare 24/iin 
excesses from young A and solar-like stars and find that the incidence is quite similar at a 
given age. Age effects would presumably cause an apparent decrease of excess incidence with 
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later type among the F, G, and K stars because F stars will evolve off the main sequence 
quickly enough to bias our sample toward younger objects (Section l6.3p . 

In the end, this discussion may still be suffering from a relatively small number of 
K stars sampled. An ongoing Spitzer program (PID 30490) to survey nearby stars that 
were not observed in other programs — a sample that includes ~400 K stars — should help 
unravel these statistics. Nonetheless, given our result, it is unlikely that future surveys will 
find a strong trend among F, G, and K stars of comparable ages — a range of spectral types 
that spans more than a factor of 50 in stellar luminosity. This behavior is counter to our 
expectations and is a challenge to models of debris disk evolution. 



6.3. Effects of age 

The lack of strong dependence on spectral type lets us combine data on various types to 
study the evolution with stellar age. Figure [T2l shows individual R24 and Xm determinations 
for the stars in our FGK sample whose ages are known, as a function of system age. There 
is no correlation apparent for these individual sources, so we look to binned data in a larger 
combined sample for evidence of trends. 



We once again take the union of the data presented here and tha t presented in iBeichman et al 



( 2006b ). but this time we include only spectral types F0-K5 from the Beichman et al. ( 2006bl ) 



sample (that is, we exclude the latest K stars). This is because the data we present in this 
paper covers the range F0-K5, and we want the best match to our combined sample. There 
are 10 stars in the F0-K5 TPF/SIM subsample that have no known ages, and 10 additional 
stars with ages less than 1 billion years. Of these 20, 2 systems have excesses (10%). Since 
this excess ratio is not significantly different from that of the overall FGK sample, and since 
the PID 30211 F stars sample is controlled for age but the PID 41 sample is not specifically 
controlled for age, we make no attempt to correct the larger sample for age. 

Figure d3] shows excess rate as a function of age for this combined sample. To zeroth 
order, there is no trend as a function of age: a constant excess rate of ~20% adequately fits 
the data, being consistent at la with the 10 Gyr data point and at 1.5a with the 8 Gyr data 
point. Furthermore, the 10 Gyr bin has only 7 targets in it, and the 8 Gyr bin has only 
33 targets in it (still a relatively small number). 

On the other hand, we note that the data shown in Figure [13] is suggestive of an excess 
rate that decreases with time through at least 8 Gyr. In this scenario, the 10 Gyr bin would 
be highly anomalous, although we note that this bin is clearly affected by small number 
statistics (2 excesses out of 7 stars). In other words, there may be a real trend and a real 
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evolution of planetary systems and debris disks even on the billion- year timescale. However, 
this may again be a manifestation of the (potential) observational bias shown in Figure 
as follows. 

The fraction of stars in a given age bin that are K stars increases for the later age bins. 
If K stars truly have fewer excesses (or fewer detectable excesses, according to observational 
biases) than other spectral types, Figures [131 and [T41 might indeed be showing a real decrease 
with increasing age, but caused not by a long-timescale evolution of planetary systems but 
by the increasing dominance of excess-deficient K stars at the oldest ages. The data present 
in our larger sample cannot distinguish between the competing possibilities of K stars pref- 
erentially lacking (detectable) disks, or of old stars increasingly lacking disks. It also remains 
to be seen whether the high excess rate for the oldest bin in Figure d3] is anything more than 
a small number statistics anomaly. 

The overall high rate of excess incidence in our samples (17% ) indicates tha t the 400 Myr 



decay timescale the drives the evolution of A star debris disks ( jSu et al.ll2006l ) cannot drive 
the evolution of debris disks in our >1 Gyr Sun-like sample. Instead, Sun-like stars appear 
to have a relatively constant incidence of 15%-20% that is not strongly dependent on age, 
but may be weakly dependent on age through a very long timescale decrease. 



6.4. Debris disks around Sun-like stars 



Because there is no difference in excess rate between F and G stars, we can combine 
our sample of F0-G9 stars into a single p opulation of "Sun-like" stars. To this sample of 
169 stars we add the 56 F and G stars from iBeichman et al.l (l2006b| ) to create an even larger 
sample of 225 Sun-like stars (213 stars at 24 microns). We refer to this merged sample of 
213 and 225 stars as the "Sun-like supersample." The excess rates for this supersample are 
4.21^5% at 24 microns and 16. 4+29% a ^ microns (Table HJ). With this large supersample, 
we can now state the debris disk incidence rate for Sun-like stars with quite good confidence 
(small error bars). 



6.5. Implications for planetary system formation 

The majority of the debris disk systems that we present here have excesses at 70 microns 
only, suggesting temperatures <100 K and therefore an inner edge to the disks. These dusty 
debris disks are likely produced by collisions within a swarm of planetesimals akin to the 
asteroid belt or Kuiper Belt in our Solar System. We interpret the cool temperatures we 
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derive for the dust in these disks as evidence of inner disk holes where the surface density of 
dust is much smaller than in the planetesimal ring, and potentially zero. This architecture is 
strongly reminiscent of our Solar System, where planets sculpt the edges of the planetesimal 
and dust belts. It may be that many of the systems we discuss here similarly have planets 
sculpting their dust distributions. In several cases, the known planets may indeed be the 
ones sculpting the inner edges of the disks (Tables [5] and [6]). Additionally, Lawler et al. (in 
prep.) have found, using IRS spectra, that the incidence of detectable levels of warm dust 
may be higher than the 4.2% we find here, indicating that dust in a region analogous to our 
Solar System's asteroid belt may also be somewhat common. 



Bryden et al.l (120071 ) show that the properties of disks around planet-bearing stars are 



unlikely to be similar to the properties of disks around stars without known planets. Briefly, 
the detection rates between these two populations are similar, but the planet-bearing stars 
are generally farther away and in more confused regions of the sky. They conclude that disks 
around planet-bearing stars are dustier (i.e., more massive), and probably more common, 
than disks around stars without known planets. Additionally, iTrilling et al.l (120071 ) recently 
showed that the excess rate for binary A-F stars is 9% and 40% at 24 and 70 /zm, respectively, 
significantly higher than our results for FGK stars and for Sun-like stars. We see in Figure [TT] 
that these excess rates are relatively high, compared to the large sample of single stars 
we present here. Combining these two studies, we find strong evidence that the presence 
of additional massive bodies (whether stellar or planetary companions) in stellar systems 
appears to promote higher excess rates. This effect may simply be dynamical - - more 
massive bodies means more stirring, more collisions, and more dust — but it is not clear 
that such a process could remain effective for billions of years. Further work is necessary to 
explain these results. 

Figure [TT] shows that the transition from the high excess rate around A stars to the 
more modest excess rate around Sun-like stars is gradual. This gradual decay is likely an 
age effect (consider the mean ages of the samples given in Figure [TTlh Our data show that 
excess rates for FGK stars decline slowly with stellar age, and it is not clear whether we are 
detecting a billion- year tail of debris disk evolution, a dependence on spectral type, and/or 
an observational bias. 



Young (< 1 Gyr) debris disks de c ay on 100-400 milli on year timescales (IRieke et al.ll2005 



Su et al.l 12006k iGorlova et al.l 120061 ; ISiegler et al.l 120071 ) . The presence of excesses around 
16% of Sun-like stars at billion year ages indicates that these old debris disks must be 
driven by a different evolutionary process than debris disks around those younger stars. One 
interpretation, using our Solar System as an analogy, is that the younger systems are still 
active in a Late Heavy Bombardment kind of dynamical upheaval. After a billion years, such 
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large scale processes likely have ceased in all but the most unusual systems. Debris disks are 
then produced from collisions within remaining planetesimal belts (e.g., our asteroid belt) 
that have been dynamically excited by the previous eon's dynamical stirring. However, there 
is no trend of fractional luminosity with age (Figure IT41) . although there is a lack of high 
fractional luminosity disks at old ages. 

The presence of a debris disk indicates that planetary system formation progressed 
at least to the planetesimal stage in a given system. There is no apparent dependence on 
spectral type across Sun-like stars for fractional luminosity (i.e., dust mass), dust distance, or 
even excess rate. This indicates that the processes that give rise to debris disk — planetesimal 
formation, dynamical stirring that produces collisions - - must equally be insensitive to 
stellar parameters (temperature, mass, luminosity). This may argue that planetary system 
formation is quite robust — able to occur in many different conditions. While none of the 
debris disks we observed are very similar to our own Solar System, there can be no question 
that the process of planetary system formation is quite common. 

7. Summary 

We observed nearly 200 FGK stars at 24 and 70 microns with the Spitzer Space Telescope 
to search for debris disks around nearby stars like the Sun. We identify excess emission, 
including a number of cases where the observed flux is more than 10 times brighter than the 
photospheric flux. We combine our data with results from several other sources to create a 
sample of more than 350 AFGKM stars. The incidence of debris disks in a large sample of 
Sun-like stars is 4.2l^'5% at 24 microns (213 stars) and 16.412 9% a t 70 microns (225 stars). 
We find that the excess rates for A, F, G, and K stars are essentially indistinguishable, but 
with a suggestion of decreasing excess rate toward the later spectral types; this may be an 
age effect. The lack of strong trend among FGK stars of comparable ages is surprising, 
given the factor of 50 change in stellar luminosity across this spectral range. We also find 
that any decline in debris disk activity over the 1-10 Gyr time frame examined is very slow. 
This result contrasts with the more rapid decay found previously for stars 0.01-1 Gyr. This 
contrast suggests that the behavior at younger ages is dominated by events analogous to our 
Solar System's Late Heavy Bombardment. How debris disks are maintained for the ensuing 
billions of years remains an outstanding question. 

We thank suggestions from an anonymous referee that helped us clarify our results. 
This research has made use of NASA's Astrophysics Data System (ADS abstract server); 
the SIMBAD databased, operated at CDS, Strasbourg, France; and of data products from 
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the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts 
and the Infrared Processing and Analysis Center/ California Institute of Technology, funded 
by the National Aeronautics and Space Administration and the National Science Foundation. 
This work is based in part on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA 
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1255094 issued by JPL/Caltech. 
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Table 1. Target information 



Name Spec. V d Age [M/H] Ref. 

type (mag) (pc) (Gyr) 



63 


F5 


7.15 


50. 


.48 


2.9 


-0 


.08 


M 


142 a 


GUV 


5.70 


25. 


.64 


2.88 





.04 


Ca 


166 


KOV 


6.13 


13 


.70 


5.01 


-0 


.30 


Ca 


693 


F5V 


4.89 


18. 


.89 


5.16 


-0 


.34 


Ca 


1237 a 


G6V 


6.7 


17. 


.62 


2.75 





.10 


Ca 


1539 


F5 


7.03 


46. 


.32 


2.9 





.03 


M 


1581 


F9V 


4.20 


8. 


.59 


3.02 


-0 


.10 


Ca 


3126 


F2 


6.92 


41. 


.51 


3.5 


-0 


.22 


M 


3296 


F5 


6.61 


47. 


.19 


2.5 





,01 


M 


3302 


F6V 


5.52 


36 


.22 


7.76 


-0 


.02 


M 


3651 a 


KOV 


5.80 


11 


.11 


5.89 





.26 


Hy 


3795 


G3V 


6.14 


28. 


.56 


7.24 


-0 


.73 


Ca 


3823 


G1V 


5.88 


25. 


.47 


5.50 


-0 


.35 


Ca 


3861 


F5 


6.53 


33 


.34 


3.4 





.03 


M 


4307 


G2V 


6.15 


31 


.86 


7.76 


-0 


.36 


Ca 


4628 


K2V 


5.75 


7. 


.46 


8.13 


-0 


.29 


Ca 


7570 


F8V 


4.96 


15. 


.05 


4.33 





.08 


Ca 


8070 


F2 


6.66 


48 


.69 


2.5 


-0 


.20 


M 


8574 a 


F8 


7.8 


44. 


.15 


7.24 


-0 


.23 


M 


9826 a 


F8V 


4.09 


13 


.47 


6.31 


-0. 


,03 


Ca 


10476 


K1V 


5.20 


7. 


.47 


4.57 


-0. 


.20 


Ca 


10647 a 


F9V 


5.52 


17. 


.35 


6.3 


-0. 


.11 


M 


10697 a 


G5IV 


6.29 


32 


.56 


7.41 





.10 


Ca 


10800 


G2V 


5.89 


27. 


.13 


7.37 





.03 


M 


13445 a 


K1V 


6.17 


10 


.95 


5.56 


-0. 


.21 


Ca 


13555 


F5V 


5.24 


30 


.13 


2.74 


-0 


.18 


Ca 


14412 


G5V 


6.34 


12. 


.42 


3.31 


-0 


.53 


Ca 


14802 


G2V 


5.19 


21 


.93 


6.76 


0. 


.00 


Ca 


15335 


GOV 


5.91 


30 


.79 


7.76 


-0. 


.22 


Ca 


15798 


F5V 


4.75 


25 


.82 


3.16 


-0. 


.30 


Ca 



-26- 



Table 1 — Continued 
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Spec. 


V 


d 


Age 


[M/H] 


Ref. 




type 


(mag) 


(pc) 


(Gyr) 






1 fi1 fin 


K3V 

IVO V 


89 

O. OZj 


7 91 

1 . Zj X 


K 4 b 


-fi 94 

U . Ut: 


Hv 


1 70^1 a 

1 i \JO X 


nnv 

VJU V 


^ 4fl 


1 7 94 


9 4^ 


04 
-u. u^± 




1 7Q9^ 

1 i a^O 


iv i v 


fi nn 


in 


n iq 

u. iy 


1^ 

-\J. LO 




1Q^7^ 

X UO 1 o 


cnv 


4 s ) 

t: . \JO 


lU.utJ 


^ 89 


0^ 


Ca 


1 QQQ4 a 


F8V 

I O V 


u.UU 


99 38 

Z(Z( .(JO 


3 ^ 


n 09 


f!a 


9n3fi7 a 

ZjUOU 1 


VJU 


fi 41 

U . t: X 


97 1 3 

Z/ f . 1l) 




-D m 

U . U X 


M 

1VX 


906^0 


v < -J V c 


4 8^ 


Q 16 


9Q 

U. z,c/ 


01 

- U. U X 




ZjU I uu 


(49V 


^ ^4 

O. O'-k 


19 19 

X Z, . X z. 


^ 90 

tj. ZjU 


-0 90 

U. ZjU 


C!a 




m v 

V < X V 


^ 94 


1 9 D8 


7 88 


-fl 9D 

U. ZjU 


f!a 

V_ CI 


99484 


F8V 


4 98 


1 fi 84 


8 39 

O. uz 


n nn 

u.uu 


f!a 

v_ cl 


93n79 a 


V < U V 


7 1 

1 . X 


34 fi 


7 89 


-n 19 

U. 1Z 


M 

1VX 


93 c i9fi a 


F8 

I O 


7 94 


^1 98 

O X . C/O 




n n4 

U. U1 


M 

1VX 


9fi993 


VjrUI V 


fi 33 


91 1 9 

Zj X . X u 




n 1 7 

U.l 1 


f!a 

v_ CI 




K9Tva 


4 44 


18 93 

XO .ZjO 


fi fi 


n 99 

U . — — 


f!a 


981 8^ 

ZjO XOtJ 


V.! •) 


7 81 

1 . X 


39 ^fi 


3 4^ 


n 9^ 

U. ZjtJ 


UU 


3949^ 


m v 


^ ^n 

O.OvJ 


1 3 39 


1 39 

X . uZ 


n in 

u. xu 


f!a 




FfiV 

X U V 


3 19 

O. X U 


8 03 


1 

X . OO 


n ni 

U. U X 


M 

1VX 


332fi2 

-.7 -.7 Zj V z* 


F7V 

X l V 


4 72 

t: . I Zj 


1 1 65 

X ± . \JO 


3 ^2 


-n 23 

W . ZjO 


Ga 


33 ^fU 3 - 


FfiV 

I U V 


1U 


9fl 98 

ZjU . C/O 


3 48 


-n 1 1 

W . X X 


M 

1VX 


33fi3fi a 


VTU 


7 Dfi 


98 fi9 


3 94 


-n 14 

U . -Lt: 


M 


O /I A 1 1 

o4411 


Lil.olV- V 


4. f U 


1 O P. K 

12.00 


V/^ 

0. f 


n no 
-U.U8 


oa 


34721 


GOV 


5.96 


24.93 


6.17 


-0.25 


Ca 


35296 


F8Ve 


5.00 


14.66 


3.76 


0.00 


Ca 


37394 


KlVe 


6.23 


12.09 


0.49 


-0.20 


Ca 


39091 a 


G1V 


5.67 


18.21 


5.609 


0.04 


M 


40979 a 


F8 


6.75 


148.37 




-0.03 


M 


43162 


G5V 


6.37 


16.69 


0.37 


-0.16 


Hy 


43834 


G6V 


5.09 


10.15 


7.61 


0.00 


Ca 


50554 a 


F8V 


6.86 


31.03 


4.68 


-0.12 


M 


50692 


GOV 


5.76 


17.27 


4.47 


-0.11 


M 
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Table 1 — Continued 



Name Spec. V d Age [M/H] Ref. 

type (mag) (pc) (Gyr) 



52265 a 


GO 


6.30 


28.07 


6.03 


0. 


.21 


Ca 


52711 


G4V 


5.93 


19.09 


4.77 


-0. 


16 


Ca 


55575 


GOV 


5.55 


16.86 


4.57 


-0. 


.28 


Ca 


57703 


F2 


6.78 


44.09 


2.3 


-0. 


.07 


M 


58855 


F6V 


5.37 


19.90 


3.6 


-0. 


31 


Ca 


60912 


F5 


6.91 


46.43 


2.9 


-0. 


.04 


M 


62613 


G8V 


6.56 


17.04 


3.09 


-0. 


20 


Hy 


63333 


F5 


7.11 


43.07 


5.5 


-0. 


39 


Ca 


68456 


F5V 


4.76 


21.39 


2.43 


-0. 


.36 


Ca 


69830 a 


KOV 


5.95 


12.58 


4.68 


-0. 


03 


Ca 


69897 


F6V 


5.10 


15.67 


3.55 


-0. 


.26 


Ca 


70843 


F5 


7.04 


46.53 


1.6 


-0. 


01 


M 


71148 


G5V 


6.30 


21.79 


4.68 


0. 


.07 


M 


71640 


F5 


7.41 


44.92 


5.7 


-0. 


09 


M 


72905 


G1.5V 


5.65 


13.85 


0.42 


-0. 


.08 


Ca 


75616 


F5 


6.99 


35.61 


4.8 


-0. 


.17 


M 


75732 a 


G8V 


5.95 


12.53 


6.46 


0. 


.40 


Ca 


76151 


G3V 


6.00 


11.29 


1.84 


-0. 


02 


Ca 


77967 


FO 


6.60 


40.7 


3.5 


-0. 


37 


M 


79392 


F2 


6.760 


50.53 


1.8 


-0. 


.24 


M 


80218 


F5 


6.629 


39.20 


6.4 


-0. 


.28 


Ca 


82943 a 


GO 


6.54 


27.46 


4.07 


0. 


32 


Ca 


83451 


F5 


7.14 


48.47 


3.8 


-0. 


.14 


M 


83525 


F5 


6.91 


48.38 


5.3 


-0. 


.06 


M 


84117 


GOV 


4.94 


14.88 


4.24 


-0. 


.14 


M 


84737 


G0.5Va 


5.10 


18.43 


11.75 


0. 


.04 


Ca 


86147 


F5 


6.72 


45.11 


2.8 


0. 


03 


M 


88230 


K2Ve 


6.61 


4.69 


4.68 


-0. 


93 


Ce 


88984 


F5 


7.30 


51.39 


4.4 


-0. 


29 


M 


90839 


F8V 


4.83 


12.82 


3.39 


-0. 


.18 


Ca 
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Name 


Spec. 


V 


d 


Age 


[M/H] 


Ref. 




type 


(mag) 


(pc) 


(Gyr) 






9S081 


F5 

X O 


7 1 1 

1 .11 


51 3fi 


3 5 


-0 93 

U. ZjO 


M 




r u v 


5 94 

tj .Zv^i: 


O X . O^t 


O. xvj 


07 


Hp 


cfO xZjO 


vxU V 


^ in 


1 3 Q1 




01 


v_.'CL 


QQ1 26 

<J <J X ZjVJ 


F5 

1 tj 


VJ . <JVJ 


47 37 


5 7 

*J. 1 


-0 16 

VJ. XVJ 


M 

1VX 


1 ooofi7 

X UUUU I 


F5 

X O 


7 30 


40 77 


5 4 


-0 34 

VJ. ul 


M 

1VX 


1 01 959 

X VJ X ZiO o 


Gfi/8V 

VJU / O V 


fi 49 

VJ . t:Zj 


fi4 79 

Ut:. I Zj 


1 9fi 

XVJ. C/U 


-0 fiO 

vj. vjvj 


VF 

V X 


1 ni 501 

X VJ X tJVJ X 


G8Ve 


5 39 

tj .OZ- 


9 54 


1 1 9 

X . X Zj 


03 

VJ .VJO 


v_ ; d 


1 02438 


G5V 


6 48 


1 7 77 

X 1 . 1 I 




-0 3fi 

VJ . OVJ 




1 09870 

X vjzjO I VJ 


F9V 

X U V 


3 61 

O.U1 


1 90 

X VJ. C/VJ 


4 47 


90 

VJ .ZjVJ 


V_ (I 


X WO I I o 


F5 

X <j 


fi 87 

U.O 1 


49 95 


9 4 

Zj . t: 


Ofi 

VJ .VJVJ 


M 

1VX 


1 04731 

X Wt: 1 O X 


F6V 

X VJ V 


5 15 

*J . X tj 


94 90 


1 83 

X . (JU 


-0 21 

VJ . Zj X 


Ca 

V_ (1 


1 0498 5 a 


G9TTT 


5 80 

*J . (J W 


109 04 

X VJZj . VJt: 




-0 18 

VJ . X(J 


M 

1VX 


1 0591 2 

X \JxJ U X _ 


F5 

X tj 


6 9 s ) 


50 95 

tJVJ . ZjtJ 


1 8 

X . O 


-0 04 

VJ . VJt: 


M 

1VX 


1 0975fi 


F5 
X o 


fi 98 

VJ . C/O 


4fi 55 

t:U. tJtJ 


3 4 

■ J . L ± 


-0 99 

VJ . ZjZj 


M 

1VX 


1 10897 


gov 


fi 00 

U .vjvj 


1 7 37 

X I . O 1 


9 7 


-0 59 


Ga 




G7V 

VI ( V 


fi 31 


17 17 

X I . X I 


1 23 

X . Zj*J 


1 8 

VJ . X (J 


Hv 


1 1 1 


F^ 


fi 9 


48 45 


1 7 


-VJ. Uz, 


M 


1 1 21 64 

XX— IUt: 


G1 V 

VI X V 


Fi 89 

<J . O <J 


39 73 


3 43 

*J . TI- J 


24 

VJ . Zjt: 


Ca 


1 14fi1 3 

X XtiVJ X O 


G3V 

VI - J v 


4 85 


90 48 


5 97 

U. Zj 1 


Ifi 

VJ . X VJ 


VF 

V X 


1 1471 n 


F9 5V 


4 9fi 

t: . Zj VJ 


9 15 

<J . XtJ 


2 29 

Zj . Zj <J 


Ofi 

VJ . VJVJ 


Ca 

V_ (A 


11/1 *7ona 

114729 


G3 V 


6.69 


35.00 


6.76 


-0.26 


r V 


114783 a 


K0V 


7.57 


20.43 


4.37 


-0.11 


Hy 


115383 


GOV 


5.22 


17.95 


0.40 


0.04 


Ca 


115617 


G5V 


4.74 


8.53 


6.31 


-0.03 


Ca 


117043 


G6 


6.50 


21.34 




0.22 


Hy 


117176 a 


G5V 


5.00 


18.11 


5.37 


-0.09 


Ca 


118972 


Kl 


6.93 


15.60 




-0.03 


Hy 


120005 


F5 


6.50 


44.92 


2.9 


0.05 


M 


120136 a 


F7V 


4.50 


13.51 


1.91 


0.30 


Ca 


120690 


G5V 


6.45 


19.92 


2.24 


-0.11 


Ca 
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Name 


Spec. 


V 


d 


Age 


[M/H] 


Ref. 




type 


(mag) 


(pc) 


(Gyr) 






1 99869 

X ZjZjO WZj 


vjtzj .ox V 


6 09 

W . WZj 


98 68 

Z-O .UO 


6 1 1 

U.ll 


-fl 1 1 

U . X X 


M 

1VX 


1 936Q1 


F9 

1^ Zj 


6 8fl 
u.ou 


OU .Oct 


1 ^ 

X . o 


fl HQ 


M 

1VX 


l^UUUU 


F7V 

X ( V 


4 i n 

t: . X W 


14 ^7 

XH:.0 I 


Zj . 1 W 




C!a 


1 97334 

X Zj 1 -X UT: 


G5V 

VI O V 


6 40 

W .t:W 


93 ^7 

ZjO . O 1 


6 99 

W . <_/Zj 


s ) 

w . wo 


Ca 


IZiOO X X 


IVU 


7 ^1 

1 .O X 


16 ^7 

1U .O I 




n D8 

w. wo 


Hv 


i srufin 

1JUt:UU 


F^i 

X o 


7 93 


48 31 

4:0 .O X 


3 ^ 

O. O 


-fl D3 

w . wo 


M 

1VX 


1 3flQ48 


VI Z. V 


^ 88 


1 7 Q4 


n 87 


n 90 

W. ZjW 




1 33009 

X l)UU WZj 


F9V 

X V 


Fi 64 


43 33 

TT-X . OO 


9 4^9 

Zj . t:0 U 






1 34fl83 


X O V 


4 Q3 

H:. C/O 


1Q 79 

XC/ . 1 Zj 


1 73 

X . I O 


n nn 

W. WW 


C!a 


1 34Q87 a 


gw 

VI -J v 


6 4^ 

W .t:0 


9^ 6^ 

ZjO .UO 


7 76 


n 36 

w. ow 


v_ a 


1 36064 

X O WW Wt: 


F8V 

X O V 


5 10 


95 31 

ZjO . O X 


4 64 

t: . Wt: 


-0 10 

W . X w 


Ca 


1361 18 a 


F8 

X (J 


6 94 


52 27 

UZ . Zj 1 




-0 15 

W . X xj 


M 


1 41 1 28 

14:1 X ZjO 


F5 

X O 


7 01 

1 . W X 


51 02 

O X . WZj 


O 


-0 94 

W . Zjt: 


M 

1VX 


149373 

XH:ZjO 1 O 


FQV 

X U V 


4 69 

t: . WZj 


xo .oo 


8 1 3 

O . X o 


-fl 40 


C!a 


1 4986f) 

Xt:ZjO WW 


F6TV 

X Ul V 


3 8^ 


1119 

X X . X Zj 


9 88 

Zj .OO 


-fl 1 3 

W . XO 


C!a 

V_ (A 


i 43i o^ 


F5 

X O 


6 76 

W . 1 W 


46 08 


3 3 

o . o 


-0 10 

W . X w 


M 

1VX 


1 43761 a 

14:0 1 U X 


VxW V 


^ 40 

O .t:W 


1 7 43 

X 1 . t: O 


7 41 

I .t: X 


-fl 96 

W .ZjW 


C!a 


XH:OW 1 O 


KDV 


6 67 

U.U 1 


1 8 1 ^ 

XO . XO 


6 Q9 

U . C/Zj 


fl ^(] 

W. OW 


C!a 


146933 


m v 

VI X V 




1 ^ 36 

XO .ou 


4 ^7 

t:.0 1 


fi n ^ 

w. WO 


f!a 

V_ (A 


14Q661 


K9V 

1V_ V 


^ 76 


Q 78 

£7. / O 


1 1 7 


fi nn 


Hv 


1 room 

15zo91 


G8 V 


6.64 


16.94 


0.58 


-0.18 


Hy 


154088 


G8IV-V 


6.59 


18.08 


5.89 


0.30 


Hy 


157214 


G2V 


5.40 


13.57 


6.46 


-0.41 


Ca 


160691 a 


G3IV-V 


5.15 


15.28 


6.67 


0.16 


Ca 


166620 


K2V 


6.37 


11.10 


5.01 


0.07 


Hy 


168151 


F5V 


5.03 


23.50 


2.53 


-0.17 


Ca 


168443 a 


G5V 


6.92 


37.88 


8.51 


0.10 


Go 


169830 a 


F8 


5.91 


36.32 


7.24 


0.13 


Ca 


171886 


F5 


7.18 


49.65 


3.7 


-0.33 


M 


173667 


F6V 


4.20 


19.09 


3.39 


-0.01 


Ca 
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Name Spec. V d Age [M/H] Ref. 

type (mag) (pc) (Gyr) 



176441 


F5 


7.08 


44. 


.72 


3.7 


-0. 


.27 


M 


177830 a 


K0 


7.18 


59 


.03 


8.5 


0. 


36 


Go 


181321 


G5V 


6.49 


20 


.86 


0.5 








181655 


G8V 


6.31 


25. 


.23 


4.57 


0. 


02 


VF 


185144 


KOV 


4.70 


5. 


61 


3.24 


-0. 


.14 


Hy 


186408 


G1.5V 


5.96 


21 


.62 


10.4 


0. 


.08 


Ca 


186427 a 


G3V 


6.20 


21 


.41 


7.41 


0. 


.08 


Ca 


188376 


G5V: 


4.70 


23 


.79 


4.47 


-0. 


.13 


Ca 


189567 


G2V 


6.07 


17. 


.71 




-0. 


30 


Ca 


190007 


K4V 


7.48 


13 


.11 


1.5 b 








190248 


G7IV 


3.56 


6. 


.11 




0. 


32 


Ca 


191408 


K3V 


5.31 


6. 


.05 


7.88 


-0. 


.58 


Ca 


193664 


G3V 


5.93 


17 


.57 


4.7 


-0. 


.18 


M 


196050 a 


G3V 


7.6 


46. 


.93 




0. 


.23 


FV 


196378 


F8V 


5.12 


24. 


.20 


6.02 


-0. 


30 


Ca 


196761 


G8V 


6.37 


14. 


.65 


4.27 


-0. 


60 


Hy 


197692 


F5V 


4.15 


14. 


.67 


1.94 


0. 


.04 


Ca 


200433 


F5 


6.92 


48. 


.36 


1.4 


-0. 


.09 


M 


202884 


F5 


7.27 


42. 


.05 


5.2 


-0. 


.29 


M 


203608 


F8V 


4.22 


9. 


22 


10.196 


-0. 


.64 


Ca 


206860 


GOV 


6.00 


18. 


.39 


5.00 


-0. 


.20 


Ca 


207129 


GOV 


5.58 


15. 


.64 


5.76 


-0. 


.15 


Ca 


209100 


K4.5V 


4.69 


3. 


63 


1.38 


0. 


.04 


Ca 


210277 a 


GOV 


6.63 


21 


.29 


6.76 


0. 


.22 


Ca 


210302 


F6V 


4.92 


18. 


.74 


5.37 


0. 


.06 


Ca 


210918 


G5V 


6.26 


22 


.13 


3.85 


-0. 


.18 


Ca 


212330 


G3IV 


5.31 


20 


.49 


7.94 


0. 


00 


Ca 


212695 


F5 


6.95 


51 


.10 


2.3 


-0. 


.05 


M 


213240 a 


GOV 


6.80 


40. 


.75 


2.77 


0. 


.04 


M 


216345 


G8III 


10.21 


100 
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Name Spec. V d Age [M/H] Ref. 

type (mag) (pc) (Gyr) 



216437 a 


G2.5IV 


6.06 


26.52 


7.15 


0.10 


Ca 


216803 


K4V 


6.48 


7.64 


0.2 b 


0.07 


Sa 


217014 a 


G4V 


5.49 


17.12 


7.41 


0.05 


Ca 


217813 


G5 


6.66 


24.27 


0.71 


-0.07 


M 


219134 


K3V 


5.56 


6.53 




0.20 


Hy 


219983 


F2 


6.67 


47.71 


4.7 


-0.27 


M 


220182 


Kl 


7.36 


21.92 


0.26 


0.01 


Hy 


221420 


G2V 


5.82 


31.76 


5.52 


0.37 


So 


222143 


G5 


6.58 


23.12 




0.08 


Hy 


222368 


F7V 


4.13 


13.79 


3.93 


-0.18 


Ca 


222404 a 


K1IV 


3.23 


13.79 




-0.05 


Ca 


222582 a 


G3V 


7.70 


41.95 


6.63 


0.02 


Go 


225239 


G2V 


6.10 


36.79 




-0.50 


Ca 



''Star with known planet (s). 



b Age from iBarnea (120071 ). 



Note. - Spectral types and V magnitudes are from B06 
and from SIMBAD an d Hipp arcos. Ages are from B06 and 
from Nordstrom et al.l (120041 ) except as indicated. To date, 



Barnes! ( 2007 ) has derived ages for ~15% of the stars in this 
sample, and in general the agreement with our compiled liter- 
ature values is good. When the number examined becomes a 
significant fraction of the total sample, we will re-examine 
excesses as a function of age using the new, systematic 
age definitions (future work). Distances are from Hipparcos 
(jPerryman et al.lll997l ). A number of stars have no good ages, 
and a handful have no good metallicities. Metallicities are 
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(VF) 



(FV) Fischer fc Valentil (120051) : (Go) iGonzelaz et al 



from the indi cated references: ( Ca) ICavrel de Strobel et al 

fhooi ) 

(2001) 
Jl995 1 



(Sa) ISantos et al 



fHv) iHavwoodl (20011): (M) Marsakov fc Shevelev 



J200l|); (So) 



Sousa et al. 



(j200fih : 



Valenti fc Fischerl (120051 ). 
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Table 2. Observing log 



Name 


Int. time 


Int. time 


AORKey a 


PID 




24 /im 


70 /im 








(sec) 


(sec) 






Do 


A S 
4o 


000 


1 7Q/I RQn/l 
1 ( 0400U4 


Qnoi 1 
OUZl l 


1/19 
14Z 


n 
u 


1 nan 
iuyu 


1 971 c;nns 

1/ 1 lOUUo 


A 1 
41 


1 49 

14Z 


4o 


OOO 


4n891 7fi 
4UoZl 1 O 


41 
41 


1 P.P. 

100 


/I8 
4o 


a /in 

44U 


/i n^n79n 

4UOU ( ZU 


/1 1 
41 


0»O 


/IS 
4o 


/i /in 

44U 


a n f ?nQ7fi 
4uouy i o 


A 1 
41 


1 9^7 
I/O ( 


/IS 
4o 


1 ^nn 

10UU 


a nRnfi79 
4UOUO ( z 


/1 1 
41 


100y 


/IS 
4o 


9Q1 
ZOl 


1 7Q/I 

1 ( 040ZOZ 


^noi 1 

OUZl 1 


1 ^81 
lOol 


/IS 
4o 


1 9fi 
1ZO 


A n^l 9^9 
4U01ZOZ 


A 1 
41 


^1 9fi 
01Z0 


/IS 
4o 


9Q1 
ZOl 


1 7 f ? f ?7finn 
1 ( OO ( ouu 


^noi 1 

OUZl 1 


ozyo 


/IS 
4o 


1 Q9 
10Z 


1 7Q/infi79 
1 1 04UO ( Z 


^noi 1 

OUZl 1 


OOUZ 


/IS 
4o 


/i/in 

44U 


/i n^i /i ss 

4U014oo 


A 1 
41 


0001 


/IS 
4o 


9Q1 
ZOl 


/in^i 7/1/1 

4U01 1 44 


/1 1 
41 


o ( yo 


/IS 
4o 


oou 


a n^onnn 

4UOZUUU 


A 1 
41 


OoZO 


/I S 
4o 


/i /in 

44U 


/in f ?99c;fi 

4UOZZOO 


/1 1 
41 


OoOl 


A S 
4o 


1 n7 


1 ( oooooz 


^noi 1 

OUZl 1 


A ^07 
40U ( 


/I S 
4o 


oou 


a nQ9t;i 9 

4UOZ01Z 


A 1 
41 


40Zo 


AS 
4o 


1 9fi 
1ZO 


4n^97fi8 
4UOZ 1 Oo 


41 
41 


( ( u 


/IS 
4o 


9Q1 
ZOl 


/in'i'ino/i 

4UOOUZ4 


/1 1 
41 


sn7n 

oU ( U 


/IS 
4o 


9Q1 
ZOl 


1 7Q/1 1 fiQfi 

i i 04ioyo 


^noi 1 

OUZl 1 


8574 





1007 


4005376 


41 


8574 


92 


2181 


8775936 


41 


9826 


48 


231 


4033280 


41 


10476 


48 


126 


4033536 


41 


10647 


136 


336 


7865856 


717 


10697 


48 


881 


4033792 


41 


10800 


48 


440 


4034048 


41 


13445 


48 


126 


4034304 


41 


13555 


48 


440 


4034560 


41 


14412 


48 


881 


4034816 


41 



-34- 



Table 2 — Continued 



Name 


Int. time 


Int. time 


AORKey a 


PID 




24 /iin 


70 fim 








(sec) 


(sec) 






1 48D9 


48 


1 9fi 


4035079 


41 

'-t X 




48 


fi50 


4035398 

'-tUOOO^O 


41 

t: X 


1 5798 

XtJ 1 UO 


48 


440 


4035584 


41 

t: X 


1 fi1 fifl 


48 


931 


4035840 


41 

t: X 


X i WtJ X 


48 


355 


4036096 


41 

t: X 


1 7Q95 


48 


1 9fi 


403R359 


41 

t: X 


1 Q373 


48 
'-to 


1 9fi 


4:UtJUUUO 


41 

t: X 




48 


44fl 


4D8nfi4n 


41 

t: X 


9(13fi7 


48 

'-tO 


i 300 


1971 1Q3fi 

XZj I X X jOU 


41 

t: X 




48 


1 9fi 


403fi8fi4 


41 

t: X 


9f)7fifi 


48 


931 


40371 9(1 

4:1)0 I X ZjU 


41 

t: X 


90807 


48 

'-tO 


44fl 


403737fi 


41 

t: X 


99484 


48 


931 


4(1(184(1(1 


41 

t: X 


9307Q 


48 
'-to 


91 81 

X O X 


4080384 


41 

t: X 


935QR 


48 


91 81 


1971 1 1 fi8 

IZi 1 XXX uo 


41 

t: X 


9RQ93 


48 


i 300 


40fiOQ98 


41 

t: X 


27442 


48 


1 26 


4081 152 


41 

t: X 


981 8 5 


1 3fi 


91 81 

X O X 


4005fi39 


41 

t: X 


30495 


48 


231 


4037632 


41 


oUooz 


/I o 

48 


zol 


/I 097000 

4Uo / 888 


/I 1 

41 


33262 


48 


126 


4038144 


41 


33564 


48 


231 


4038400 


41 


33636 





1300 


4059136 


41 


33636 


48 


2181 


8777728 


41 


34411 


48 


126 


4038656 


41 


34721 


48 


650 


4038912 


41 


35296 


48 


231 


4061184 


41 


37394 


48 


126 


4039168 


41 


39091 


48 


355 


4039424 


41 
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Name 


Int. time 


Int. time 


AORKey a 


PID 




24 /iin 


70 jum 








(sec) 


(sec) 






4DQ7Q 


48 


21 81 


1 971 1 (S8D 

-LZ, t _L 1UOU 


41 


43162 


48 


1 007 


4039680 


41 


43834 


48 


1 26 


4039936 


41 




n 

u 


i 3no 


408(11 28 


41 

t: X 


OUOOt: 


48 


21 81 


8775680 


41 

t: X 


50fiQ9 


48 


44D 


4(14(11 Q2 


41 

t: X 


^226^ 


n 


1 DQD 


1271 5776 

XZj i XO I I u 


41 

t: X 


^226^ 


48 


1 DQD 


4D62Q76 


41 

t: X 


^271 1 

OZi I X X 


48 


44D 


404(1448 

t:Ut:Ut:t:0 


41 

t: X 


OOO I o 


48 


355 

OOO 


4(14(17(14 

1Ut:U 1 Ut: 


41 

t: X 


577(13 

O I I WO 


48 


231 


1 733888(1 

X I OuOOOU 


3(121 1 

OW^< X X 


58855 
ooooo 


48 


OOO 


4n4DQ6(l 


41 

t: X 


6DQ1 2 


48 


231 


1 7342Q76 


3021 1 

OVJ^t X X 


6961 3 


48 


1 f)Qf) 


4(141 21 6 

t:Ot: IZi X \J 


41 

t: X 


63333 

\JOOOO 


48 


231 


1 73381 12 

X I OOO X XZj 


3021 1 

OW^< X X 


68456 


48 


231 


4D41472 

t:Ot: 11 ( Z( 


41 

t: X 


69830 


48 


355 


4041 728 


41 


6Q8Q7 


48 


231 


4(1(16144 

tiWWVJ X X X 


41 

t: X 


7(1843 


48 


231 


i 7344nnn 


3(121 1 


71 1 /I O 

f 1148 


/I o 

48 


881 


AC\A 1 HO /I 

41141984 


A 1 

41 


71640 


48 


355 


17346560 


30211 


72905 


48 


355 


4042240 


41 


75616 


48 


231 


17336064 


30211 


75732 


48 


440 


4042496 


41 


76151 


48 


545 


4042752 


41 


77967 


48 


126 


17336832 


30211 


79392 


48 


231 


17342720 


30211 


80218 


48 


107 


17336576 


30211 


82943 


48 


1636 


4063232 


41 
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Name 


Int. time 


Int. time 


AORKey a 


PID 




24 /iin 


70 /iin 








(sec) 


(sec) 






834^1 


48 


931 


1 73494R4 


3091 1 


83^9^ 


48 

40 


931 


1 7341 nc >9 

X 1 04:1 UO 


3091 1 

tJU^ X X 


841 1 7 

(Jt:1 X 1 


48 


196 


4043008 

4W4*J uu(j 


41 

4 X 


84737 


48 

40 


931 


40439R4 


41 

4 X 


861 47 


48 

40 


1 39 


1 734041 fi 

X 1 Ot:Ut:1U 


3091 1 

tJyJZi X X 


88930 


48 

40 


931 


404377fi 

4U40 1 1 U 


41 

4 X 


88Q84 


48 

40 


3^ 


1 7347840 

X 1 04 ( Ot:U 


3091 1 

tJU^ X X 


Q083Q 


48 

40 


931 


40066"i6 


41 

4 X 


Q3081 


48 

40 


3^ 


1 734^094 


3091 1 

tjyjZi x x 


Q4388 


48 

40 


440 

111) 


4044039 

4L/44L/0 £j 


41 

4 X 


Qol 98 


48 

40 


931 


4044988 


41 

4 X 


QQ1 9fi 


48 

40 


931 


1 7341 440 


3091 1 

<0\JZj X X 


X UUUU I 


48 

40 


3^ 


1 7337088 

X I OO I woo 


3091 1 

<0\JZj X X 


1 01 9^Q 


48 
40 


1 0Q0 


4044^44 

4:1/4:4:04:4: 


41 

4 X 


1 ni 

X U XtJU X 


48 

40 


931 


4044800 

4U44QUU 


41 

4 X 


1 09438 


48 

40 


1 0Q0 


404 c >0 I i6 


41 

4 X 


1 0987(1 


48 

4U 


196 


404^)31 9 

4U4(JU1ai 


41 

4 X 


1 03773 

X KJO I I O 


48 

40 


931 


1 7340Q98 


3091 1 

tJU^ X X 


1 04731 


48 


440 

1 1 V 1 


404^894 


41 


1 n /i no e 
104985 


48 


a A n 
440 


1 071 01 no 

1/ i Izlyz 


A 1 

41 


105912 


48 


355 


17345280 


30211 


109756 


48 


231 


17342208 


30211 


110897 


48 


755 


4046080 


41 


111395 


48 


755 


4046336 


41 


111545 


48 


355 


17344768 


30211 


112164 


48 


440 


4046592 


41 


114613 


48 


355 


4046848 


41 


114710 


48 


126 


4047104 


41 


114729 


48 


1510 


4084480 


41 
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Name 


Int. time 


Int. time 


A0RKey a 


PID 




24 /iin 


70 /iin 








(sec) 


(sec) 






1 1 478^ 




9181 




41 

^± X 


X A-OOOO 


48 


931 

Zv*_> X 


4047360 


41 

t: X 


X X tJ\J X 1 


n 


1 321 


12712960 


41 

t: X 


X XtJU X l 


48 


1 9fi 

X Zv\J 


404761 fi 


41 

t: X 


X X I 1)4:0 


48 


1 DD7 


4047879 

t:Ut: I O 1 Zv 


41 

t: X 


I I 71 7fi 

III 1 1 U 


n 

u 


44D 

ill) 


19716800 

X Zv i 1UOUU 


41 

t: X 


1 1 71 7fi 

XXI 1 1 u 


48 


1 9fi 

X Zv\J 


40481 98 

lUiO XZvO 


41 

t: X 


1 18Q79 

X XOC7 1 Zv 


48 


44n 

ill) 


4048384 


41 

t: X 


1 9000^ 


48 


1 07 


1 733Q648 


3091 1 

tJUZ X X 


1 9D1 ^fi 

X ZvW x ovj 


48 


1 9fi 

X Zv\J 


4048640 


41 

t: X 




48 


1 DQD 


40488Q6 


41 

t: X 


1 998fi9 


48 


^4^ 


4D4Q1 ^9 


41 

t: X 


1 9^fiQ1 

X ZjkJXJU X 


48 




1 7344^1 9 

X i u 1 1(J X Zv 


3091 1 

O WZv X X 


X ZiUUUU 


48 


1 9fi 

X Zv\J 


4006Q1 9 

IUUUjI Zv 


41 

t: X 


1 97^S4 

X Zv I OOI 


48 


1 DD7 

XUU I 


4D4Q4fl8 


41 

t: X 


1 98^1 1 

X ZiOu X X 


n 
u 


1 DQD 


40834^6 


41 

t: X 


1 2831 1 

X ZiUu X X 


48 


2181 


87761 92 

U 1 I \J X U Zv 


41 

t: X 


X OU4:UU 


48 




1 7347^84 


3091 1 

tJUZ X X 


1 ^f)Q48 


48 


44fl 

1 1 V 1 


404Q664 


41 


1 oonno 
looUUz 


4o 


a A n 
44 U 


a a a noon 

4U4yyzu 


/I 1 

41 


134083 


48 


231 


4050176 


41 


134987 


48 


1090 


4050432 


41 


136064 


48 


231 


4050688 


41 


136118 


92 


2181 


12712448 


41 


141128 


48 


355 


17345792 


30211 


142373 


48 


126 


4050944 


41 


142860 


48 


126 


4051200 


41 


143105 


48 


132 


17340160 


30211 


143761 


48 


355 


4051456 


41 
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Name 


Int. time 


Int. time 


AORKey a 


PID 




24 /iin 


70 /iin 








(sec) 


(sec) 






i^oyj t o 


48 


1 DD7 


4051 71 9 


41 

t: X 


1 46933 


48 

tO 


9^1 

ZjO X 


4081 664 


41 

t: X 


1 49661 


48 

TiO 


196 


4051 968 

tzWO X JUU 


41 

T: X 


X O^jOU X 


48 

tiO 


1 son 


4D61 44D 

TlW\J X 1 11) 


41 

t: X 


1 ^4088 


48 

tO 


i soo 


4D61 6Q6 


41 

t: X 


1 ^7914 

X O 1 Zill 


48 

tO 


9^1 

Z<0 X 


40071 68 

1UU I xuo 


41 

t: X 


1 606Q1 
x uuuy x 


48 

tiO 


1 96 


4061 Q59 

t:UU1 JiJ^i 


41 

t: X 


X UUUZiU 


48 

tiO 


9^1 

ZjO X 


4059994 


41 

t: X 


1 681 ^1 

X UO XO X 


48 

tiO 


OOO 


4059480 

TlUOZiTlOU 


41 

T: X 


1 (S844S 

X VJ (J 11 O 


48 

tiO 


1 96 


4059736 

tiUOZj I ou 


41 

t: X 




48 

tO 


^4^ 

OtiO 


4063488 


41 

t: X 


1 71 886 

X 1 X oou 


48 

tiO 


3^ 

ooo 


1 7346048 

X 1 04:1)114:0 


3091 1 

OVJ^i X X 


1 7^667 


48 

tiO 


1 96 


4059QQ9 


41 

t: X 


1 7R441 

X f V J 11 X 


48 

tiO 


3^ 

ooo 


1 733Q3Q9 

X i OOzjOz) 


3091 1 

OVJ^i X X 


1 778S0 


48 

tiO 


9181 


4053948 

4UtJOZi40 


41 

t: X 


1 81 391 

X O X 0£t X 


48 

tiO 


1 son 


4069908 


41 

t: X 


1 81 6^5 

X (J X \JO O 


48 


755 

i OO 


4053504 

t:W OOO Wt: 


41 

T: X 


1 8^144 

X OO X tit: 


48 

tiO 


1 96 


4053760 

4UJO 1 UU 


41 

t: X 


1 86408 


48 


755 

t OO 


405401 6 


41 


18642 r 


45 


7F rr 
(00 


4U04U16 


A 1 

41 


188376 


48 


126 


4054272 


41 


189567 


48 


545 


4054528 


41 


190007 


48 


1007 


4082432 


41 


190248 


48 


126 


4054784 


41 


191408 


48 


126 


4055040 


41 


193664 


48 


440 


4055296 


41 


196050 


48 


2181 


12711424 


41 


196378 


48 


231 


4055808 


41 


196761 


48 


881 


4056064 


41 
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Name 


Int. 


time 


Int. 


time 


AORKey a 


PID 




24 


/xm 


70 


/im 








(sec) 


(sec) 






197692 




48 




126 


4056320 


41 


200433 




48 
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a Further details of each observation, including 
pointing and time and date of observation, can be 
queried from the Spitzer Data Archive at the Spitzer 
Science Center using the specified AORKey. 

Note. - HD 10647 was observed in PID 717, a 
MIPS in-orbit check-out (IOC) program. 
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Table 3. Photospheric predictions and photometry for all sources 
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Table 3 — Continued 
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Table 3 — Continued 
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Table 3 — Continued 



Name 


V 


K 


F24 


P24 


R24 


F70 
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Table 3 — Continued 
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P24 
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Table 3 — Continued 



Name 


V 


K 


F24 


P24 R24 
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(mag) 
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-0.89 



Beichman et al.l (j2006al ) found that this system has a weak excess in the 30-34 micron 



region that is not evident in our 24 or 70 micron data. 



McGregor] (J199J) 



Johnson et al.l (119681 ) 



d Derived from B-I 



Alonso et all (119941 ) 



E Derived from B-V 



^Glassl (119751 ) 



feouchet et all (|l99lh 



Johnson et al.l (119661 ) 
^Veederl Jl974h 
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Persson et al.l (119771 ) 



Alonso et all (Il998f ) 



' ICarney fc Aaronsonl (119791 ) 



qGlassI fll974h 



c Mould & Hvland 


(1976) 


'Haas & Leinert 


c 


L990) 



Note. 



V mag nitudes are from Hipparcos, with typical errors 0.01 mag 



(IPerrvman et al. 



( jSkrutskie et al. 



1997). K magnitudes are from 2MASS, with typical errors 0.02 mag 



2006h . except as indicated. All of these sources are strongly detected at 



24 microns, with intrinsically large S/N. We calculate P70 as F24/9.10, except for the 
six cases with 24 micron excesses and P24 (see text), where we calculate P70=P24/9.10. 
For HD 207129, a seventh star with 24 micron excess, we identify the excess through 
Kurucz fitting; see text. The implied excess for HD 191408 is spurious (see text). There 
are 12 stars for which K band data is not available, and therefore no P24 or R24 deter- 
minations. However, even for these stars, we can still find P70 as F24/9.10, and measure 
X7o- For these twelve stars, X70 is actually a lower limit, since there could be small 
24 micron excesses that we cannot measure. 



Table 4. Summary of excess rates 



Sample 


24 micron 


24 micron 


70 micron 


70 micron 




Number 


excess rate 


Number 


excess rate 






(%) 




(%) 


A 


30 


7 (+3/-3) 


27 


26 (+10/-7) 


FGK sample 


184 


3.8 (+1.7/-1.2) 


196 


16.3 (+2.9/-2.8) 


FGK sample, planets 


45 


6.7 (+5.8/-2.1) 


48 


20.8 (+7.0/-4.6) 


FGK sample, no planets 


139 


2.9 (+2.2/-0.8) 


148 


14.9 (+3.4/-2.5) 


All F 


110 


7.3 (+3.3/-1.2) 


117 


17.9 (+4.1/-3.0) 


All G 


103 


1.0 (+2.2/-1.0) 


108 


14.8 (+4.0/-2.8) 


All K 


51 


5.9 (+5.2/-1.8) 


51 


13.7 (+6.2/-3.5) 


M 


62 


0.0 (+2.9) 


13 


0.0 (+12) 


Sun-like supersample (FG) 


213 


4.2 (+2.0/-1.1) 


225 


16.4 (+2.8/-2.9) 



Note. - The A stars sample includes stars with ages >600 Myr from ISu et al.l (120061 ). 
The FGK sample is all data presented in this paper, from PIDs 41 and 30211, plus the 
three additional stars listed in Section 15.11 The "All F /G/K" samples are pro duced by 
combining data presented in this paper with data from iBeichman et al.l (1 2006b'). for th e 
given spectral types. The M stars results are taken directly from iGautier et al.l (120071 ). 
The Sun-l ike supersample i s defin ed in Section 16.41 as all F and all G stars from this paper 
and from IBeichman et al.l (j2006bl ). combined. In all of these samples, all stars that have 
excesses at 24 microns also have excesses at 70 microns (except some of the A stars, where 
only upper limits are available at 70 microns for some sources). In all cases in this paper, 
we cite binomial errors that i nclude 68% of th e prob ability (equivalent to the la range for 
gaussian errors) , as defined in iBurgasser et al.l (120031 ) . 
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Table 5. Color temperature for excesses detected at both bands 



Name 


F24 

L ^^excess 


F70 

J- 1 ^excess 


T color 




fd 




(mJy) 


(mJy) 


(K) 


(AU) 


xl0~ 5 


166 


20 


89 


87 


9.1 


5.9 


3126 


6 


104 


66 


21.8 


13 


10647 a 


33 


842 


62 


21 


34 


69830 b 


74 


9 


1900 


0.02 


560 








400 


1.0 


20 


101259 


30 


8 


271 


0.98 


6.0 


105912 


16 


29 


110 


7.7 


7.9 


207129 c 


24 


274 


72 


15.3 


12 



a This star has a known planet orbiting at 2.03 AU 
fautler et al.ll2006h . 

b The first line for HD 69830 here shows our naive solu- 
tion for this system, based on the fluxes at 24 and 70 mi- 
crons, and assuming that all of the flux arises from contin- 
uum emission. A dust grain temperature of nearly 2000 K 
is not sensible, though, as this is above the sublimation 
temp erature for any reasona ble dust composition. How- 



ever. 



Beichman et al.l (j2005bl ) found that a significant frac- 



tion of the excess flux from this system appears in emission 
features, meaning that our continuum assumption overes- 
timates the dust temperature (with consequences for dust 
distance and fractional luminosity). The second line for 
HD 6 9830 here shows the bes t valu es for these param eters 
from iBeichman et al.l (j2005bl ) and iLisse et al.l (120071 ) , de- 
rived from more appropriate model fitting of the emission 
spectrum. Note that this star ha s three planets orbiting 
within 0.65 AU Lovis et al.ll2006h . 



c There is no good K magnitude available for this star, 
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but Kurucz model fitting clearly indicates an excess at 
24 microns. The excess at 70 microns is indicated by 
comparison to the Kurucz model and also by the ratio 
F24/F70. 



Note. - T dust is found by fitting a black body to the 
measured excesses. This black body is then scaled to fit the 
observed fluxes, and the ratio of this scaled black body's 
integrated flux to the star's integrated flux gives the frac- 
tional luminosity (f d ). We solve for R dust using Equa- 
tion [TJ e Eri is not listed here, but does have excess at 
both b ands (Backman et al., in prep.) and a planet orbit- 
ing it jliatzes et alihoool ). 
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Table 6. Dust properties for systems with excess at 70 microns only 



Name 


Teff 
(K) 


Tdust 
(K) 


Rdust 
(AU) 


Max f d 
xlO^ 5 


Min/ d 
xlO^ 5 


1581 


6100 


218 


2.3 


1.6 


0.2 


3296 


6440 


89 


17 


3.8 


2.5 


17925 


5080 


103 


4.8 


4.4 


2.2 


19994 a 


6200 


168 


4.1 


3.5 


0.7 


20807 


5950 


150 


4.3 


1.5 


0.4 


22484 


6200 


142 


5.8 


4.3 


1.2 


30495 


5950 


88 


12 


3.0 


2.0 


33262 


6300 


158 


5.0 


1.1 


0.2 


33636 b 


6030 


81 


16 


3.5 


2.6 


50554 c 


6200 


79 


19 


3.8 


2.9 


52265 d 


6030 


92 


12 


4.0 


2.4 


57703 


6890 


78 


22 


5.0 


3.8 


72905 dc 


5900 


123 


6.2 


1.5 


0.6 


t OV) 1 u 




77 


94- 


^ 4 


9 (S 


76151 f 


5830 


120 


5.9 


1.0 


0.4 


82943 g 


6030 


69 


22 


10.0 


8.8 


110897 


6030 


94 


12 


2.3 


1.4 


115617 


5770 


97 


8.3 


3.3 


1.9 


117176 11 


5770 


153 


3.4 


7.7 


1.8 


118972 


5080 


101 


5.0 


4.6 


2.4 


128311 j k 


5250 


106 


5.1 


2.7 


1.3 


206860 1 


6030 


124 


6.6 


1.5 


0.6 


212695 


6440 


75 


24 


6.4 


5.1 



(Mayor et al. 


2004) 


b This star has a 


(Butler et al. 


2006) 
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c This star has a known planet orbiting at 2.28 AU 
fautler et alfaod ). 

d This star has a known planet orbiting at 0.5 AU 
fautler et allbood ). 



Beichman et al.l (j2006al ) derive a fractional luminosity 



of 2.9 x 10 5 based on 70 micron photometry and IRS 
spectroscopy. 



Beichman et al.l (j2006al ) derive a fractional luminosity 



of 1.2 x 10 5 based on 70 micron photometry and IRS 
spectroscopy. 

g This s tar has two known planets orbiting at 0.752 and 
1.19 AU jMavor et al.lbooi iLee et al-lbood ). 

h This star has a known planet orbiting at 0.484 AU 
(feutler et al1l2006h . 



Beichman et al.l (j2006al ) found no excess for this star 



in IRS spectra out to 35 microns, but calculate a frac- 
tional luminosity of 2.6 x 10~ 5 based on 70 micron pho- 
tometry (the same data that are presented here). 

k This s tar has two know n planets orbiting at 1.10 and 
1.76 AU JVoet et al.lboOol ). 



Beichman et al.l (j2006al ) derive a fractional luminosity 



of 0.7 x 10 5 based on 70 micron photometry and IRS 
spectroscopy. 



Note. — These dust temperatures are maximum tem- 
peratures, and the distances consequently are minimum 
distances. Both maximum and minimum fractional lu- 
minosity (fd) values are given; the former corresponds to 
the dust temperatures and radii given in this table (see 
text), while the latter corresponds to 50 K dust tem- 
perature solutions. These two values give the range of 
fractional luminosities that fit the data (see Figure [8]). 
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Fig. 1. — Distribution of spectral types for the 193 stars presented here. Hashes from lower 
right to upper left show stars from PID 30211, and hashes from lower left to upper right 
show stars from PID 41. Most of the targets in PID 30211 are F5 stars (difficult to see in 
this representation). 
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Fig. 2. — Histogram of distances for the 193 stars presented here. Three stars from PID 41 
have distances off this plot: two around 100 pc, and one around 150 pc (Tabled]). As in 
Figure [U hashes from lower left to upper right show stars from PID 41, and hashes from 
lower right to upper left show stars from PID 30211. 
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Fig. 3. — Histogram of ages for the 193 stars presented here. As in previous figures, hashes 
from lower left to upper right show stars from PID 41, and hashes from lower right to upper 
left show stars from PID 30211. 
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Metallicity 



Fig. 4. — Metallicities for the 189 stars presented here with known metallicities (see Tabled]). 
As in previous figures, hashes from lower left to upper right show stars from PID 41, and 
hashes from lower right to upper left show stars from PID 30211. 
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Fig. 5. — Histogram of R24 values, as described in the text. The shaded histogram shows 
the total number of stars in each bin, while the dotted region shows the subset for which 
K magnitudes are not derived from 2MASS. There is no systematic effect by using K mag- 
nitudes from other sources. The vertical line at 1.10 shows our excess threshold, as defined 
in the text. The core of this distribution is symmetric about unity. Seven systems have R24 
values greater than 1.10; one of these is spurious, but the others are real excesses. See text 
for discussion. 
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Fig. 6. — Distribution of X70 values. The same data are presented in both panels. The top 
panel shows a smaller range in X70 in order to show the shape of the core of the distribution. 
The lower panel shows data for all 193 systems. Vertical lines in each panel show X70 — 3.0, 
our threshold criterion. The bin sizes in the two panels are not the same (0.5 and 1 in top 
and bottom, respectively). Two noteworthy systems are labeled in the lower panel. The 
two "no K" systems HD 19994 { X7 o = 4.25) and HD 207129 (x 7 o = 24.6) are included in 
this figure. HD 100067 has a very low X70 value of -4.43, a result of a negative F70 value 
for this star, which means that this star was not detected. No stars that were detected have 
Xro < -3.0. 
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Fig. 7.— R24 and %70 for the 181 stars with P24 (Table [3]). The R24 and X70 excess 
thresholds are shown by vertical and horizontal lines, respectively. The box in the lower left 
shows the area enclosed by (-3a,3a) in the two dimensions. Targets within this box show the 
core of our sample, having no excess at either band. In addition to showing the distribution 
of our sources in this parameter space, we also confirm here that no obvious systematic 
effects derive from our technique for deriving P24 and P70 exist. Such a systematic error 
would be manifest as a correlation between R24 and X70 i n the (core) non-excess sample 
(if the photospheric models were incorrect, for instance). Instead, we see that this core 
distribution is essentially a random scatter distribution (fully uncorrelated) , confirming that 
our technique for predicting the 24 and 70 micron fluxes is appropriate. The regions with 
excesses at 24 only, 70 only, and at both bands are labeled. Targets with excesses at both 
bands are labeled (HD 10647 is off the plot to the top). HD 191408 clearly shows no excess at 
70 microns, and we find that its 24 micron excess is spurious (see text). No other 24 micron 
excesses are spurious, since all of the rest have clear 70 micron excesses as well (excepting 
HD 69830, though the excess around that star is clear for many other reasons; see text). We 
also label the other five systems with the largest %70 values, e Eri (not observed as part of 
the programs presented here; see Backman et al. (in prep.)) is also plotted and labeled; it 
is clear that neither the 24 or 70 micron excesses for this famous debris disk are remarkably 
large. The star's proximity (3.2 pc) accounts for its disk's brightness and prominence. 
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Fig. 8. — Derived fractional luminosities and dust distances as a function of stellar effective 
temperature (i.e., spectral type). We show maximum (downward triangles) and minimum 
(upward triangles) fractional luminosities in the top panel, and minimum dust distances 
(upward triangles) in the lower panel. For the seven systems with excesses at both wave- 
lengths we can solve for the fractional luminosity and dust distance exactly (not limits), 
so these values are shown as diamonds (green symbols). These seven systems have rela- 
tively high fract i onal lu mino sities. For HP 69 830, we plot here the published values from 
Beichman et al.l (j2005bl ) and iLisse et al.l (120071 ). As a point of comparison, our Solar Sys- 
tem's asteroid belt resides at 3-5 AU, and our Kuiper Belt exists at 30-50 AU. The lack of 
excesses at large distances and small fractional luminosities around the latest stellar types 
may be an observational bias. 
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Fig. 9. — Range of possible fractional luminosities as a function of dust temperature for 
three noteworthy systems. The top right of each panel (hot, bright emission) is generally 
ruled out by our 24 micron measurements. Submillimeter observations could place limits on 
the top left of each panel (bright, cold emission). Based on the measured 70 micron excess, 
possible dust temperatures and luminosities are shown as shaded regions (gray for 3cr error 
limits, black for la limits). The approximate characteristics of the asteroid and Kuiper belts 
are shown for comparison. The plausible range o f dust propert i es for HD 693 overlaps with 
the properties of the dust in our Kuiper Belt (e.g.. lTeplitz et al.l (Il999l )). The dust properties 
for HD 3126 are well-constrained due to its detection at both 24 and 70 microns. 
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Fig. 10. — Spectral energy distribution (SED) for HD 207129. The total SED is shown by 
the dashed purple curve; the components are stellar (black dashed line) and dust (solid black 
line) emission. There is no good K magnitude for this star (the 2MASS K band error is large), 
but we can find a satisfactory Kurucz fit, as shown. Visible (Hipparcos) and near infrared 
(2MASS) photometry are shown as blue circles and purple diamonds, respectively. IRAS 
measurements are shown as red symbols. Our MIPS photometry is shown as blue diamonds. 
We find an excess at both 24 and 70 microns and therefore derive the temperature and 
fractional luminosity given in the legend. The high IRAS 25 micron flux is due to another 
bright source ~30 arcsec away. 
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Fig. 11. — Exc ess rate a s a fun ction of spectral type (for "old" s tars, >600 Myr). Th e A stars 
values are from lSu et al.l (120061 ) and the binaries values are from lTrilling et al.l (120071 ) (a study 
of debris disks in A3-F8 binaries; different shading is used here to emphasize that the sample 
is defined differently, i.e., includes multiplicity as a requirement ) . The F , G, and K samples 
are the union of the data presented here and in iBeichman et al. (h006bl ). The mean ages for 
each of the A, F, G, and K samples are given at the top of the figure (each value has a large 
standard deviation of perhaps 2-3 Gyr; additionally, ages are notoriously hard to derive for 
main sequence FGK stars). While formally these data are consistent with no dependence 
on spectral type, a decrease in excess as a function of spectral type is weakly suggested. 
However, this could easily be an age effect (see text and Figure [T2"]) . The excess rate for 
(presumably old) M sta rs is 0%, with upper limits (binomial errors) of 2.9% at 24 microns 
and 12% at 70 microns JGautier et al.ll2007h (Table 0}. 
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Fig. 12. — Individual R24 and X70 determinations as a function of age for all stars in our 
FGK sample with good ages (164 have good ages and K band photometry and therefore 
R24 measurements; 175 have good ages and X70 measurements). No trend is apparent. The 
solid horizontal line in the top panel indicates R24=1.0, where there is no excess. (The 
no excess case in the bottom panel would be X70 — 0-0, which cannot be displayed on this 
logarithmic scale.) The horizontal dashed lines indicate our thresholds for excess (R24=1.10 
and X70 = 3.0). 
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Fig. 13. — Excess rate as a function of age for the F0-K5 sample (our data plus that from 



Beichman et al.l (j2006bl )). The horizontal error bars show not the uncertainties in stellar 
age (which could be large) but rather the bin widths. The data could be consistent with 
essentially no trend as a function of age (that is, with one point low and another point 
low), but are suggestive of an overall decrease with age. This could indicate that debris 
disks continue to evolve even on billion-year timescales. However, this trend could also be a 
manifestation of the stellar components of these bins, since the older bins will be increasingly 
dominated by K stars (see Figure [TTj) . The number of stars in the six bins are as follows, 
from youngest to oldest: 24, 57, 60, 52, 33, 7. Since the oldest bin has only 7 stars in it 
(2 with excesses), the high value may be a small number statistical anomaly. 
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Fig. 14. — Fractional luminosi ty as a function of age for the excesses in the "all FGK" 
sample (our data plus that from iBeichman et al.l (j2006bl )). The individual data points show 
fractional luminosity determinations from Tables [5] and [6] (using the maximum fractional 
luminosity values for the 70 micron excesss only cases). The larger grey symbols show the 
means in the same age bins as used in Figure [T3l The horizontal error bars show the bin 
widths, and the vertical error bars show the la errors on the means. There is no obvious 
trend of fractional luminosity with age, though there is a lack of high fractional luminosity 
disks older than 6 Gyr. The systems with the three largest fractional luminosities have been 
omitted: e Eri (0.3 Gyr, 29), HD 69830 (4.7 Gyr, 20), and HD 10647 (6.3 Gyr, 34) (where 
the fractional luminosities are given in units of 10 -5 ). Including these three obvious bright 
outliers does not change the result that no trend is evident, and just produces much larger 
error bars on the mean values. Our Sun has an age of 4.5 Gyr, and our Solar System's 
fractional luminosity is 10 6 10 7 . 



